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Multijet cross sections at the LHC and Tevatron are sensitive to several distinct kinematic energy 
scales. When measuring the dijet invariant mass rrijj between two signal jets produced in association 
with other jets or weak bosons, mjj will typically be much smaller than the total partonic center- 
of-mass energy Q, but larger than the individual jet masses m, such that there can be a hierarchy 
of scales m <^ rrijj <^ Q. This situation arises in many new-physics analyses at the LHC, where 
the invariant mass between jets is used to gain access to the masses of new-physics particles in a 
decay chain. At present, the logarithms arising from such a hierarchy of kinematic scales can only 
be summed at the leading-logarithmic level provided by parton-shower programs. We construct an 
effective field theory, SCET+, which is an extension of soft-coUinear effective theory that applies 
to this situation of hierarchical jets. It allows for a rigorous separation of different scales in a 
multiscale soft function and for a systematic resummation of logarithms of both rrijj /Q and m/Q. 
As an explicit example, we consider the invariant mass spectrum of the two closest jets in e^e~ 3 
jets. We also give the generalization to pp ^ N jets plus leptons relevant for the LHC. 



I. INTRODUCTION 

A detailed understanding of events with several 
hadronic jets in the final state is of central importance at 
the Large Hadron Collider (LHC) and the Tevatron. This 
is because the Standard Model and most of its exten- 
sions produce energetic partons in the underlying short- 
distance processes, which appear as collimated jets of 
energetic hadrons in the detectors. 

Multijet processes are inherently multiscale problems. 
For generic high-pr jets, there are at least three rele- 
vant energy scales present: the pr or total energy of 
a jet, which is of order the partonic center-of-mass en- 
ergy Q of the collision, the invariant mass m of a jet, 
which is typically much smaller than Q, and Aqcd, the 
scale of nonperturbative physics in the strong interaction. 
To separately treat the physics at these different energy 
scales one has to rely on factorization theorems. With a 
factorization theorem in hand, the long-distance physics 
sensitive to Aqcd can be determined from experimental 
data, while short-distance contributions can be obtained 
through perturbative calculations. 

For most multijet events, the jets will not be equally 
separated and they will not have equal energy. Instead 
there will be some hierarchy in the invariant masses be- 
tween jets or the jet energies due to the soft and collinear 
enhancement of emissions in QCD. Depending on the 
measurement, the cuts imposed to enhance the sensitivity 
to physics beyond the Standard Model introduce sensi- 
tivity to these additional kinematic scales. For example, 
requiring a large pt for the leading signal jet(s) and a 
smaller px for additional jets introduces sensitivity to a 
new kinematic scale, namely the pr of the subleading jet. 
Another example is the measurement of the dijet invari- 
ant mass of jets, which is used to identify new particles 
decaying to jets, a special case being the identification 
of boosted heavy particles by measuring the invariant 



mass of two subjets within a larger jet. As a final ex- 
ample, consider the CDF measurement + 2 jets [ij, 
which shows an excess in the invariant mass of the two 
jets around rrijj ^ 150 GeV. The excess is dominated by 
back-to-back jets, where the subleading jet is rather soft 
with Pt ~ 40 GeV, whereas the total invariant mass of 
the Wjj system is of order Q ^ 300 GeV. The challenges 
in this type of measurement are clear: the recent DO mea- 
surement 0] finds no evidence for an excess. A precise 
theoretical understanding of multijet processes with mul- 
tiple scales is clearly valuable. 

Such configurations with multiple different kinematic 
scales present are not necessarily well described in fixed- 
order perturbation theory. The sensitivity to different 
scales can give rise to Sudakov double logarithms of the 
ratio of those scales, such as {rrijj /Q) or \n^{pT /rrijj), 
at each order in perturbation theory. If the scales are 
widely separated, the logarithmic terms become the dom- 
inant contribution and the convergence of fixed-order 
perturbation theory deteriorates. A well-behaved per- 
turbative expansion that allows for reliable predictions is 
obtained by performing a resummation of the logarithmic 
terms to all orders in perturbation theory. At present, 
the resummation of such kinematic Sudakov logarithms 
in exclusive multijet processes can only be carried out at 
the leading-logarithmic level using parton-shower Monte 
Carlos. 

In this paper, we consider exclusive jet production in 
the kinematical situation where all jets have compara- 
ble energy ^ Q and two of the jets become close^, as 



We fondly refer to this kinematic region as "ninja", since it cor- 
responds to Ui ■ rij <^ 1, where and Uj are the light-cone 
directions for nearby jets. We thank Iain Stewart for publiciz- 
ing this codename so that this kinematic region has now become 
known as the ninja region. 



2 




(a) All jets equally separated. (b) Two jets close to each other. 

FIG. 1: Different kinematic situations and relevant scales for the case of three jets with invariant mass m. On the left, the 
invariant masses between any two jets are comparable, Sij = 2qi ■ qj ~ Q^, and the only relevant scales are Q, m, and m} jQ. 
On the right, the dijet invariant mass between jets 1 and 2, t = 2q\ ■ q2, is parametrically smaller than that between any other 
pair of jets, so there are two more relevant scales, y/i and rr? j^ft. 



illustrated in Fig. 1(b) for the case of three jets. We are 
interested in the dijet invariant mass mjj between the 
two close jets, which is much smaller than the other di- 
jet invariant masses of order Q, but much larger than 
the invariant mass m of the individual jets, i.e., there 
is a hierarchy of scales m ^ mjj ^ Q. In this case, 
the cross section contains two types of logarithms, those 
related to the mass of the jets, ln^(m/Q), as well as kine- 
matic logarithms ln^(mjj/(3). For m^j ^ Q, all jets are 
well separated, as in Fig. l(a)[ and the jet-mass loga- 
rithms ln^(m/Q) in the exclusive jet cross section can 
be resummed [S-Q using soft-collinear effective theory 
(SCET) HO]- 

In this paper, we construct a new effective theory, 
SCET_|_, which is valid in the limit m <^ rrijj Q. The 
added complication in this case arises from the fact that 
one needs to separate the soft radiation within a given jet 
from the radiation between the two close jets, giving rise 
to two different scales. In regular SCET, both of these 
processes are described by the same soft function, which 
therefore contains multiple scales. Soft functions with 
multiple scales have been observed in SCET before, and 
it has been suggested that this requires one to "refac- 
torize" the soft function into more fundamental pieces 
depending on only a single scale. This was first pointed 
out in Ref. [5|]. Here we explicitly construct for the first 
time an effective theory that accomplishes a refactoriza- 
tion of the soft sector and separates different scales in a 
soft function. Using SCET+, we derive the factorization 
of multijet processes in the limit m <^ nijj <C Q, where 
each function in the factorization theorem depends only 
on a single scale. The renormalization group evolution in 
SCET+ then allows us to sum all large logarithms aris- 
ing from this scale hierarchy, including those in the soft 
sector. 

It is worthwhile to note that the multijet events we con- 
sider in this paper are part of a broader class of kinematic 
configurations that give rise to multiple disparate scales. 
The case we address here of small dijet invariant masses 
belongs to the class of configurations for which the kine- 



matics of the final-state jets introduces additional kine- 
matic scales. In our case this gives rise to large logarithms 
of ratios of dijet masses \T\{mjj /Q). Other configurations 
which give rise to large kinematic logarithms, such as 
those with a hierarchy of jet prs, may require a different 
effective-theory treatment, which we leave to future work. 
These kinematic logarithms are in contrast to so-called 
"nonglobal" observables [12], which introduce additional 
scales by imposing parametrically different cuts in differ- 
ent phase space regions. This corresponds for example 
to a hierarchy between individual jet masses mi <C rrij^ 
giving rise to logarithms of the form \D.(mi/mj). The 
structure of such logarithms has been recently explored 
using SCET in Refs. 0,111. 

In the next section, we explain the physical picture 
of the effective-theory setup. In Sec. IIIIl we discuss the 
construction of SCET-|_ , which requires a new mode with 
coUinear-soft scaling to properly describe the soft radia- 
tion between the two close jets. As an explicit example of 
the application of SCET+, we consider the simplest case 
of e^e~ 3 jets, for which in Sec. IIVI we derive the fac- 
torized cross section in the limit m <C mjj <^ Q, and in 
Sec.|V]we obtain all ingredients at next-to- leading order 
(NLO). In Sec. |Vl we also discuss the consistency of the 
factorized result in SCET_(-, and show how the usual 3-jet 
hard and soft functions in SCET are separately factor- 
ized into two pieces each. Readers not interested in the 
technical details of this example can skip over Sees. IIVI 
and |Vl In Sec. I VII we generalize our results to the case 
oi pp N jets plus leptons. In Sec. I VIII we present nu- 
merical results for the dijet invariant mass spectrum for 
e^e~ — > 3 jets with all logarithms of m/Q and nijj /Q re- 
summed at next-to-leading logarithmic (NLL) order. We 
conclude in Sec. I Villi 
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II. OVERVIEW OF THE EFFECTIVE FIELD 
THEORY SETUP 

Effective fiefd tlieories provide a natural way to sys- 
tematically resuni large logarithms of ratios of scales ap- 
pearing in perturbation theory. This is achieved by inte- 
grating out the relevant degrees of freedom at each scale. 
The renormalization group evolution within the effective 
theory is then used to sum the logarithms between the 
different scales. SCET provides the appropriate effective- 
theory framework to resum the logarithms arising from 
coUinear and soft radiation in QCD. 

For later convenience, for each jet we define a massless 
reference momentum 

g,f ^Ei{l,ni) = \qiT\{coshr]i,ni±,sinhr]i) , (2.1) 

where Ei and fii are the energy and direction of the ith 
jet, and (jiT, rji, and j_ = qiTl\<iiT \ are its transverse mo- 
mentum, pseudorapidity, and transverse direction. The 
four- vector (l,rij) is called n^. 

Given two light-cone vectors rii and rij with nf — n'j — 
0, we decompose any four-vector into light-cone compo- 
nents 



— rij ■ p ■ 



+ ni- p- 



Ui ■ rij 



■Pa 



(2.2) 



The subscripts on J-ij specify with respect to which light- 
cone directions the perpendicular components are de- 
fined. To simplify the notation we will mostly drop them, 
unless there are potential ambiguities. To discuss the 
scaling of a four-vector with respect to the rii direction, 
we use 



{ni-p,ni-p,pj_^) = {p+,p ,pj_), 



(2.3) 



where hi stands for any other light-cone vector in a di- 
rection that is considered parametrically different from 
rii, i.e., parametrically rii ■ hi 1. In terms of these 



light-cone coordinates we have p 



p+p- 



-Pa 



A. Equally Separated Jets 

We first review the situation for equally separated jets, 
as depicted for three jets in Fig. l(a)[ where we show the 
relevant energy scales. Formally, this case is defined by 
considering the pairwise invariant masses between any 
two jets to be parametrically of the same typical size Q. 
In addition, the invariant masses of all jets are paramet- 
rically of the same typical size m. Hence, we have the 
scaling 



2qi ■ Qj ^ 



(2.4) 



where Pi denotes the total four-momentum of the jets. 
For m much smaller than Q, we have the hierarchy of 
scales 



MS 



Q 



<^ Hj C:^ m <^ fiH — Q ■ 



(2.5) 



The cross section contains logarithms scaling like 
\n{m^/Q^). To resum these in the effective theory, one 
first matches QCD onto SCET at the hard scale Q; see 
Fig. 1(a) In this step, one integrates out all degrees of 
freedom with virtualities ^ Q^. The relevant modes in 
SCET (technically SCETj) below that scale are collinear 
and ultrasoft (usoft) modes. They have momentum scal- 
ing 

P,^g(A^l,A), pus-^Q{X',X\X^), A = ^. 

. ^^-^^ 

(The scaling for the collinear momentum pc here is un- 
derstood to be defined with respect to its corresponding 
jet direction.) To see this, first note that to describe the 
collinear emissions within each jet we must have one set 
of collinear modes in each jet direction. Since the typi- 
cal invariant mass of the jets is m^, the collinear modes 
must have virtuality ~ Q^X^ — m? . Direct interactions 
between two collinear modes in different directions are 
not allowed, since they would produce modes with virtu- 
ality ~ , which have already been integrated out, i.e.. 



Q(A2,1,A) + Q(1,A2,A) ^Q(1,1,A). 



(2.7) 



Hence, interactions between collinear modes can only 
happen via usoft modes, which can couple to any 
collinear mode without changing its virtuality: 



Q(A2, 1, A) + Q{X\ X^, A2) ^ Q(A2, 1, A) , 



(2., 



In the next step, one integrates out the collinear modes 
at the scale m, leaving only the soft theory with usoft 
modes of virtuality ~ Q^A^ = m'^/Q^. 

In each step one performs a power expansion in the 
ratio of the lower scale divided by the higher scale. 
In the first step m/Q = A, and in the second step 
{rn^ /Q)/m — A, so the expansion parameter is the same 
in both cases. By Lorentz invariance, in the end the ex- 
pansion is actually in the ratio of the virtualities, i.e., in 
A^. There will be no power corrections of C(A). 

Using this effective theory, one can derive a factoriza- 
tion theorem for the cross section forpp — > N jets, which 
has the schematic form 0, [^, 0, [l^, [l^ 

N 
i=l 

(2.9) 

The hard matching coefficient Cn describes the short- 
distance partonic 2 ^ N process. It arises from integrat- 
ing out the hard modes at the scale Q when matching full 
QCD to SCET. The beam functions Ba^b and jet func- 
tions Ji describe the collinear initial-state and final-state 
radiation, respectively, forming the jets around the ingo- 
ing and outgoing primary hard partons. They arise from 
integrating out the collinear modes at the intermediate 
jet scale m. The remaining matrix element in the soft 
theory yields the soft function Sn- In general, Cn is a 
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vector in color space and Sn is a matrix, while the beam 
and jet functions are diagonal in color. 

Each function in Eq. (j2.9|) explicitly depends on the 
renormalization scale /i. This dependence cancels in the 
product and convolutions of all functions on the right- 
hand side, since the cross section is /i independent. Since 
the different functions each only contain physics at a sin- 
gle energy scale, they can only contain logarithms of fi 
divided by that physical scale. In general, the logarithms 
appearing in the hard, jet, and soft function are of the 
form 



Sn 



In 



In 



2 



(2.10) 



(For the iV-jet soft function this can be seen explicitly 
from the results in Ref. [3].) Hence, since all ^ and 
all rrii ^ m as in Eq. ()2.4p , there are no large logarithms 
when evaluating each function at its own natural scale. 



m = Q : 



m 



Using the renormalization group evolution in the effec- 
tive theory, each function can then be evolved from its 
own natural scale to the common arbitrary scale /x, which 
sums the logarithms in Eq. (j2.10l) . Combining all func- 
tions evolved to ^ as in Eq. (|2.9p then sums all logarithms 
of the form ln(m^/(3^) in the cross section. 



Two Jets Close To Each Other 



In the situation depicted in Fig. 1(b) the invariant 



mass of two of the jets becomes parametrically smaller 
than all the other pairwise invariant masses between jets. 
In the following, we take the two jets that are close to 
each other to be jets 1 and 2 and use t = si2 to denote 
their invariant mass, and to denote all other dijet 
invariant masses. We then have 



m <^t — si2 Si 



(2.12) 



In principle, the factorization theorem in Eq. (j2.9p can 
still be applied in this case, since the invariant masses of 
all jets are still much smaller than any of the dijet invari- 
ant masses. However, the hard matching coefficient Cjv 
now depends on two parametrically different hard scales, 
\/t and Q, and from Eq. (|2.10|) it contains corresponding 
logarithms ln(/i^/(5^) and ln(/i^/i). This means there is 
no single hard scale that we can choose that would 
minimize all logarithms in the hard matching. In par- 
ticular, choosing puj = Q as before, there are now unre- 
summed large logarithms \n{t/Q^) in the hard matching 
coefficient. 



Similarly, the soft function now depends on two para- 
metrically different soft scales, m^/^/t and rri^/Q, con- 
taining logarithms ln(/Lt^t/m'') as well as \n{iJ.^Q^ /mj). 
Hence, there is not a single soft scale /is we can choose 
to minimize all logarithms in the soft function. Choos- 
ing fis — m^/Q as before, there are still unresummed 
large logarithms \n{t/Q^) in the soft function. In the 
soft function these naturally arise as hmi-rij. 

To be able to resum the logarithms of ln(t/Q^) we have 
to perform additional matching steps at each of the new 
intermediate scales ^/t and rr? jyft^ as shown in Fig. 1(b) 



At the scale Q we match QCD onto SCET as before, in- 
tegrating out hard modes of virtuality ^ . This ef- 
fective theory has collinear modes with virtuality t and 
corresponding soft modes, 



Pc Q(Aj,l,At) , Pn 

(2.13) 

There is one set of collinear modes for each of the jets, 
except for the two close jets 1 and 2. The latter are 
described at the hard scale Q by a single set of collinear 
modes with virtuality t in a common direction nt- Since 
the total invariant mass between the two jets is t, such nt- 
coUinear modes can freely exchange momentum between 
the two close jets without changing their virtuality. This 
matching corresponds to performing an expansion in Aj. 

In the next step, at the scale v^, we match SCET onto 
a new effective theory SCET+ , integrating out all modes 
of virtuality t. Below this scale we now have separate 
collinear modes with virtuality rri^ for each jet, including 
jets 1 and 2, 



Pe = Q(A^l,A) 



A 



(2.14) 



Note that for the well-separated jets, this matching at \/t 
will have no effect, since we do not perform a measure- 
ment in those directions that is sensitive to this scale. 
This means the virtuality of their collinear modes is sim- 
ply lowered from t to m^. On the other hand, for jets 1 
and 2 we match a single rit collinear sector in SCET onto 
two independent ni collinear and n2 collinear sectors in 
SCET+. 

As before, the collinear modes cannot directly inter- 
act with each other. Interactions between the jets are 
possible via ultrasoft modes 



-Q(A2,A2,A2), 



(2.15) 



which have virtuality Q^A^ = (m^/Q)^. In addition we 
can still have collinear modes in the nt direction which 
are soft enough to interact between jets 1 and 2 without 
changing the virtuality of the ni and rt2 collinear modes. 
These collinear-soft (csoft) modes have momentum scal- 
ing in the nt direction as 



A 

a; 



m 



(2.16) 
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We will see explicitly in Sec. IV D] that these csoft modes 
are required to correctly reproduce the IR structure of 
QCD in this limit. Note also that the csoft modes are 
not allowed to couple to other collinear modes, since they 
would give them a virtuality ~ Q'^rj^ — m? /X^ ^ m? . 

To derive the csoft scaling in Eq. (|2.16p , first note that 
in order for a soft gluon to separately interact with jet 1 
or jet 2 it must be able to resolve the difference between 
the two directions ni and n2. Hence, its angle with re- 
spect to the directions rii or rt2 has to be of order of the 
separation between these two directions, so it must have 
the scaling 



Pc 



,(A?,l,At). 



(2.17) 



At the same time, for the contribution of Pcs to the in- 
variant mass of either jet 1 or jet 2 (or equivalently the 
virtuality of a ni or n2 collinear mode) to be ^ m?, its 
plus component must have the same size as a collinear 
plus component, i.e.. 



Qcs ~ Qv 



(2.18) 



which implies the scaling in Eq. (|2.16p . Hence, the rit 
csoft modes can be thought of as the soft remnant of the 
original rit collinear modes in the SCET above the scale 
Vi. Note that the csoft degrees of freedom can resolve 
the two directions ni and n2 dynamically, such that only 
a single csoft mode is required that couples to both of 
these collinear directions. We will continue to use the 
direction rit to label this mode in the rest of the paper. 
The virtuality of the csoft modes is 



m 
T 



(2.19) 



which is the intermediate soft scale we already encoun- 
tered. Hence, the main feature of SCET+ is that its soft 
sector contains two separate degrees of freedom, a csoft 
mode with virtuality m^/t, and a usoft mode with vir- 
tuality m^/Q^. In the next step, we integrate out the 
collinear modes with virtuality m? at the scale m, leav- 
ing only the soft sector of SCET+ [denoted by soft+ in 
Fig. l(b)| consisting of csoft and usoft modes. Finally, at 
the csoft scale rr? j^/t we integrate out the csoft modes, 
which leaves only the usoft modes. 

In SCET+, the cross section for pp ^ N jets takes the 
schematic form 



(2.20) 



N 



Ba{^l)Bb{^JL)\{J^{^x) ®S+{pL)®SN-l{tx) 



Here, each piece arises as the matching coefficient from 
one of the matching steps described above, and only de- 
pends on a single scale, allowing us to resum all large log- 
arithms. In effect, the additional matching at \/t factor- 
izes the original hard matching coefficient Gn in Eq. ()2.9|) 
into two pieces, Cn-i and C+, each only depending on 



Mode 


bcalmg (+, — , -LJ 


Virtuality 


collinear 


Q(A2,1,A) 


{\Qf = 


csoft 




{r]XQf = mVt 


usoft 


Q{\\X\X^) 





TABLE I: Modes in SCET4 



a single scale Q and y/t, which enables us to resum the 
logarithms ln(i/(5^) present in the original hard match- 
ing. Similarly, the additional matching step at the scale 
TV? l\/t effectively factorizes the original multiscale soft 
function S'jv in Eq. (12. 9p into two separate pieces, 5+ and 
Sat-i, each only depending on a single soft scale, vr? j \/t 
and rr? IQ, respectively. This then enables us to sum all 
logarithms ln(t/Q^) that were present in the original soft 
function. 

Note that for the hard matching coefficient, the kine- 
matic situation of t <C has been addressed before [itI - 
[l9| using a two-step matching procedure similar to our 
matching step at the scale \/t. It was shown through 
an explicit one-loop calculation that this matching sepa- 
rates the two scales and t in the hard function from 
one another. In Appendix |^ we show that this holds 
to all orders in perturbation theory using reparametriza- 
tion invariance (RPI) ^2^1 of the effective theory. In these 
previous analyses, however, the soft sector of the theory 
below was not fully considered (or assumed to be 
that of standard SCET). In this paper we give a com- 
plete description of the effective theory below the scale 
\ft including the appropriate soft sector, allowing us to 
accomplish a similar scale separation and resummation 
of the logarithms of tjQ"^ in the soft sector of the theory. 
We stress that the unresummed logarithms in the origi- 
nal Cat and Sat are equally large, so it is essential to be 
able to resum the logarithms in the soft sector. 



III. CONSTRUCTION OF SCET+ 

In this section we construct SCET_|_ , an effective theory 
containing the usual collinear and usoft modes as well as 
the new intermediate csoft mode. The modes and their 
scaling are summarized in Table HI We will show that the 
interactions between the three different types of modes 
can all be decoupled at the level of the Lagrangian by 
appropriate field redefinitions, similar to the BPS field 
redefinition to decouple collinear and usoft modes in reg- 
ular SCET. 



A. Review of Standard SCET 

In SCET, the momentum of collinear particles in the 
n direction are separated into a large label momentum p 
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and a small residual momentum k, 



n-p- 



(3.1) 



The momentum components scale asfi-p Q, p± ^ QX, 
and k ^ Q)? . The corresponding quark and gluon fields, 
£,n.,p{x) and An.p{x)^ are multipole expanded with expan- 
sion parameter A. They have fixed label momentum, and 
particles with different label momenta are described by 
different fields. Derivatives acting on the fields pick out 
the residual momentum dependence, ^ k^, while the 
large label momentum is obtained using the label mo- 
mentum operator [l(]| 



(3.2) 



When acting on several collinear fields, "P^ returns the 
sum of the label momenta of all n-coUinear fields. 

The interactions between collinear fields can only 
change the label momentum but not the collinear direc- 
tion 71, so it is convenient to define fields with only the 
direction n fixed, 



p#0 



). (3.3) 

p#0 



The sum over p here excludes the zero-bin p = 0. This 
avoids double-counting the usoft modes, which are de- 
scribed by separate usoft quark and gluon fields. When 
calculating matrix elements, we implement this by sum- 
ming over all p and then subtracting the zero-bin contri- 
bution, which is obtained by taking the limit p — >■ [2l[ . 

The Lagrangian for a collinear quark in the n direction 
in SCET at leading order in A is well known and given 

by i 



in ■Dn + gn-Aus + ilpni_ Wn 



where the collinear covariant derivatives are 



in • £)„ 



in ■ d 



- gn- Ar, 



(3.4) 



(3.5) 



The Wilson line Wn in Eq. (|3.4|) is constructed out of 
n-collinear gluons. In momentum space, one has 



Wn{x)^ ^ exp(^-^-:^n-An{x)'^ 

-perms 



(3.6) 



where the label operator only acts inside the square 
brackets. Wn sums up arbitrary emissions of n-collinear 
gluons from an n-collinear quark or gluon, which are 0(1) 
in the power counting. 

The Lagrangian for usoft quarks and gluons is identical 
to the full QCD Lagrangian written in terms of usoft 
quark and gluon fields. It cannot contain any interactions 
with collinear modes, since the usoft fields do not have 
sufficient momentum to pair-produce collinear modes. 



Because of the multipole expansion, at leading order 
in A the only coupling to usoft gluons in the collinear 
Lagrangian, Eq. p.4|) . is through n • A^s- This coupling 
is removed by the BPS field redefinition fll|, 

A(^Hx)=Y,l{x)Anix)Yn{x), (3.7) 
where 1^ is a usoft Wilson line in the direction n, 



P exp 



1.9 



dsn-Ausix + s n) 



(3.8) 



and P denotes path ordering along the integration path. 
Since Wn{x) is localized with respect to the residual po- 
sition X, we have 



Wi'Hx)^Y,l{x)Wn{x)Yn{x) 



, exp _ 
^ \n-Vn 



(3.9) 



Therefore, using Eq. ([3T7)) in Eq. together with 

Y^{x)\vn-d^gn- A^,{x)\Yn{x)^in-d, (3.10) 
eliminates the dependence of £„ on n-A^s-, 



n-Vn 



where now 



m 



.i^(o) 



iD, 



(0)m 



in • 9 - 



-gn 

9^n± ■ 



(3.11) 



(3.12) 



Hence, after the field redefinition there are no more in- 
teractions between usoft and collinear fields at leading 
order in the power counting, and the redefined fields no 
longer transform under usoft gauge transformations. 

With more than one collinear sector, there are sepa- 
rate collinear Lagrangians for each sector, which decou- 
ple from each other and the usoft Lagrangian, Cus- The 
total Lagrangian is then given by the sum 



^SCET = ^ -C^^' -I- Cii 



(3.13) 



where the ellipses denote the terms that are of higher 
order in the power counting. 



B. SCET+ 

To construct SCET4., we follow the same logic as in 
Sec. IIIBI To be concrete, we start from SCET with two 
collinear sectors along and nt that have been decou- 
pled from the usoft sector. 



-CsCET — l^n} + + 



(3.14) 
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and where the scahng of the momenta is set by At = 
\ftlQ. (Additional colhnear sectors for other well- 
separated jets are treated identically to 713.) Since the 
three sectors are all decoupled in SCET, we can dis- 
cuss them independently. When matching onto SCET-|_, 
nothing happens for the rt3-collinear and usoft modes, 
whose momentum scaling is simply lowered from At to 
A = m/Q. 

On the other hand, as explained in Sec. Ill Bl when we 
lower the scaling from At to A, the rit-collinear modes 
are separated into 7ii-collinear and n2-collinear modes, 
which cannot interact with each other any longer, plus a 
csoft mode in the rit direction. 



-A^^+AZ. (3.15) 



In SCET+ the labels rit, ni, 712 uniquely specify whether 
we are dealing with csoft or collinear modes, so we will 
mostly drop the explicit labels "cs" and "c" . 



1. Collinear-Soft Sector 

The csoft modes are just a softer version of the rit 
collinear modes. Hence, they still have a label direction 
jit and are multipole expanded, i.e., their momentum is 
written in terms of a csoft label momentum and residual 
momentum. 



yes yes 



P'is^nfp,s^+p^^^^, (3.16) 



with momentum scaling 

nt ■ Pes ^ Qrj^ , Pcs± ^ (Q'7^)At = QrjX , 

fc^ - {QTf )\^t - QA2 , (3.17) 

and an associated csoft label momentum operator, 

KAru,P^P'csU,p- (3.18) 

The Lagrangians for csoft quarks and gluons are simply 
scaled down versions of the original rit collinear quark 
and gluon Lagrangians, Cnt ■ For example, for the csoft 
quarks 



. i/T/O) T/(0) \ T/(0)ti 77/0) ''JH fiO) 



(3.19) 



where the csoft covariant derivatives are defined exactly 
as in Eq. p.l2p but in terms of csoft gluons. Note that 
the csoft modes are decoupled from the usoft modes, as 
indicated by the superscript (0), since they are obtained 
from the decoupled rif-coUinear modes of regular SCET. 
The Vif Wilson lines are the csoft version of the Wi^^ 



in Eq. dSlD, 



■perms 



(3.20) 



where T) 



(0)m 



Just like the Wn in SCET, they sum up arbitrary emis- 
sions of rit-csoft gluons from an 7T,t-csoft quark or gluon, 
which are 0(1) in the power counting, and are required 
to ensure csoft gauge invariance. 

Similarly, the csoft gluon Lagrangian follows directly 
from the collinear gluon Lagrangian in SCET, and we 
just state the result here for completeness, 

- ^Tr{ [iP(0)-+ .g4")^ ipW^-f ] 
-f 2Tr{ct'[iPf°l., [iPS)''+ff4^, c(°)]]} 
+ iTr{[iP(°),,4?'']}'. (3.21) 

fit ■ Vntnt/'^ + P^+ int ■ 9^, c„j are ghost 
fields and a is a gauge fixing parameter. 

As we will see below, csoft gluons can still couple to 
ni and n2 collinear modes but only through their uyAn^ 
and n2-Ant components, respectively. From that point 
of view, they appear more like usoft gluons, hence the 
name csoft. However, their ni component is given by 

. "l-"t _ . "l-"t . f \ A 

ni-Ant = — - — UfAnt H — rifAnt + {ni)±^-Ant± 

^ At X Qr]'^ + 1 X QA^ + At x QrjX 

+ + QA^, 

(3.22) 

and similarly for n2- The essential feature of the csoft 
scaling is that all three terms here contribute equally and 
must be kept. This is precisely the reason why despite 
being collinear modes the csoft modes are still able to 
couple to both ni and n2 collinear modes. However, since 
the csoft modes are already multipole expanded with re- 
spect to the nt direction, ni-An^ is not an independent 
component with its own power counting, but is just a 
short-hand notation for the combination in Eq. (|3.22l) . 
The same applies to n2 ■ Aj^ . 

It is important to properly implement the usoft zero- 
bin subtraction for the csoft modes, which is necessary 
to avoid double counting the usoft region. The zero-bin 
limit of the csoft modes is defined by taking each light- 
cone component with respect to rit to have usoft scaling 
^ QX^. Hence, in the zero-bin limit, only the second 
term in Eq. p.22p contributes. 



ni ■ A„ 



^ nt ■ Ant > 



(3.23) 



while the other terms are suppressed by At. 

Note that csoft fields only couple to collinear fields 
whose direction are in the same csoft equivalence class as 
nt, as discussed above. For all other collinear fields, the 
interaction with a csoft field would increase the virtuality 
of the field such that these interactions are integrated out 
of the theory. 
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2. Collinear Sectors 



C. Operators in SCET and SCET^ 



We now turn to the rii and n2 collinear modes. To be 
specific we will use ni; the discussion is identical for n2. 
In the SCET above the scale ^/t, ni and nt belong to the 
same equivalence class. ^ This means the leading-order 
Lagrangian for ni collinear quarks directly follows from 
expanding Eq. p. lip in ry, 



where the collinear covariant derivatives, D'^n} ^ and Wil- 
son line, Wi? , are as defined in Eqs. (j3.12p and (j3.9p with 
n = n\. As anticipated, the csoft modes couple to the 
ni collinear modes via ni-ALn} ~ QA^. As in Eq. p.22p . 
all components of An^ contribute equally to this coupling. 
However, below the scale \Jt^ the n\ collinear modes know 
nothing about the rit direction, so from their point of view 
the csoft modes behave just like ordinary soft modes with 
eikonal coupling in the n\ direction. In particular, just 
as in standard SCET, we can remove the coupling be- 
tween csoft and collinear modes from the Lagrangian by 
performing a field redefinition, 

(3.25) 

4r)(^)=x(°)t(xX)(.)x(°)(.), 

where the superscript (0, 0) indicates that the collinear 
fields are decoupled from both usoft and csoft interac- 
tions. Here, Xn} is now a Wilson line in the n\ direction 
built out of (usoft-decoupled) csoft gluons. 



P exp 



15 



(3.26) 



After the csoft field redefinition for n\ and n2, there 
are no more interactions between any of the sectors. The 
above discussion is not affected by additional collinear 
sectors like rij,. The Lagrangian of SCET-). thus com- 
pletely factorizes into independent collinear, csoft, and 
usoft sectors. 



(3.27) 



1=1,2 



i>3 



^ This can be understood formally using RPI [2011 , which is a sym- 
metry of the effective theory that restores Lorentz invariance of 
the full theory that was broken by choosing a fixed direction 
for each collinear degree of freedom. One can show that rii, 
and nt can all be obtained from one another by an RPI trans- 
formation, see Ref. [2^ for a detailed discussion. 



In this section we discuss how operators in SCET+ 
are constructed from gauge-invariant building blocks. As 
an explicit example, we use e~ — ^ 3 jets with jets 1 



and 2 getting close as in Fig. 1(b) since we will use it 
in Sec. IIVI For simplicity, we assume here that jets 1, 
2, and 3 are created by an outgoing quark, gluon, and 
antiquark, respectively, such that n\ = n^, n-i = n^, 
ris = riq. The operators with the quark and antiquark 
interchanged simply follow from Hermitian conjugation. 
Note that the case where the quark and antiquark jets 
get close to each other is power suppressed, so there is 
no corresponding operator in SCET+ at leading order in 
the power counting. 

The allowed operators one can construct in SCET are 
constrained by local gauge invariance. It is well known 
that using the collinear Wilson line Wnix) one can con- 
struct gauge-invariant collinear quark and gluon fields 



Xuix) = Wl{x)in{x) , 



(3.28) 



which are local with respect to soft interactions. Hence, 
we can use them to construct local collinear gauge- 
invariant operators in SCET. 



For example, for widely separated jets as in Fig. 1(a) 
we match the matrix element for e' 
QCD onto the operator 



3 jets in full 



(3.29) 



where for simplicity we neglect the Dirac structure. 
When matching QCD onto SCET in the situation with 
two close jets as in Fig. 1(b) [ we first match onto the 



SCET operator for e+e — >■ 2 jets, 

O2 = XntXr. 



(3.30) 



describing a quark and antiquark jet in the rit and n-i 
directions. Under local usoft gauge transformations, the 
fields in different collinear sectors all transform in the 
same way, so O2 and O3 are also explicitly gauge invari- 
ant under usoft gauge transformations. 

After the field redefinition in Eq. (13. 7p . we obtain cor- 
responding redefined fields x^n'ix') and i3^'^^(x) which 
are gauge invariant under both collinear and usoft gauge 
transformations. All usoft interactions are now described 
by usoft Wilson lines explicitly appearing in the opera- 
tors, e.g.. 



"-^2 Ant ^nt "3 Arig 



(0) 

"3 ' 



(3.31) 



In SCET-|_ we can use the same definitions as in 
Eq. (13.281) to define collinear fields that are gauge in- 
variant under collinear gauge transformations. The ni.2 
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collinear fields in addition transform under csoft gauge 
transformations, Unt{x), 



(x) -> Untix) Xniix) , 



(3.32) 



As discussed in Sec. IIII Bl only the ni and n2 collinear 
fields couple to csoft gluons, thus the collinear fields do 
not transform under csoft gauge transformations. There- 
fore, to form gauge-invariant operators in SCET+ we 
have to include factors of the csoft 14,t Wilson lines. For 
example, after the usoft field redefinition in SCET, an nt 
collinear quark field in SCET is matched onto SCET+ as 



Since Vn^^^ does not transform under ni^2 collinear gauge 
transformations, the right-hand side is invariant under 
both ni^2 collinear and n„j csoft gauge transformations. 
We can think of Vn^''^ here effectively as arising from the 
csoft limit of the Wilj"*^ Wilson line inside Xnt- 

Hence, when matching SCET onto SCET+ for the sit- 



V;<"^n^)<^^(^)xi°?(^)- (3.33) 



nation in Fig. |l(b)[ the SCET operator O2 in Eq. (|33T|) 
is matched onto the SCET+ operator 



(3.34) 



where the different factors in square brackets do not in- 
teract with each other. Finally, we perform the csoft field 
redefinition in Eq. p.25p to decouple the csoft fields from 
the ni,2 collinear fields, which yields 



(3.35) 



Here, each factor in square brackets now belongs to a dif- 
ferent sector in SCET_|_, and we have shown how the ad- 
joint and fundamental color indices are contracted. Since 
the different sectors are now completely factorized, we 
will drop the superscripts (0) and (0, 0) in the following 
sections. 



D. Alternative Construction of SCET+ 



When constructing SCET_|_ in Sec. IHI Bl we started 
from the usoft-decoupled version of SCET, for which the 
csoft modes arise from the usoft-decoupled nt collinear 
sector in SCET. By simply lowering the scaling in the 
usoft sector from At to A, we have implicitly used the fact 
that to be consistent and maintain the usoft decoupling 
one has to simultaneously lower the scaling of the usoft 
subtractions for the nj-collinear sector. This is the reason 
why the csoft modes arise as the csoft limit of the nt 
collinear modes. The advantage of this approach is that 



the matching onto SCET_|_ really only happens within 
one collinear sector of SCET. 

Alternatively, we can also be completely agnostic 
about the theory above the scale ^/t, and simply write 
down the Lagrangians for all the modes in SCET_|_ 
and use appropriate field redefinitions to decouple them. 
In the end, the matching calculation will ensure that 
SCET+ reproduces the correct UV physics, while having 
the right degrees of freedom ensures that the IR physics 
of the theory above is reproduced. The latter can be 
checked explicitly by testing whether the IR divergences 
in the theory above are reproduced in SCET_|_. 

This procedure should of course give the same final 
result. Since it is instructive to see how it does, we will 
briefly go through it here. The discussion for the n2 
collinear sector is again identical to that for rii, so we 
will ignore it. We start by writing down the Lagrangians 
for the collinear and csoft modes, 



Srii 



(3.36) 



where the power counting in the multipole expansion re- 
stricts the possible interactions. Note that we have added 
collinear, c, and csoft, cs, labels here for clarity. For sim- 
plicity, we only write down the quark Lagrangians and 
drop the perpendicular pieces, indicated by the ellipses, 
which are not relevant for this discussion. Note that the 
csoft gluons, A'^^ , only couple to the ni (and 712) collinear 
modes, while the usoft gluons, Aus, couple to all collinear 
sectors as well as the csoft sector. 

We first perform the usual usoft field redefinition in 
Eq. (|3.7p on all three sectors. 



A^^^Hx) = Ylix)A^MY,,M- 



(3.37) 



As far as the coupling to usoft gluons is concerned, the 
csoft sector is just another collinear sector, so we get 

£(0) _ Fc{0) \- . r)c(O) , yt y .4cs(0)yt y 

J lcc(O) 

r Sni ' 



^"3 



tc(0) 



(0) 



(3.38) 



where both the csoft and the 713 collinear sectors are now 
decoupled from the usoft. 

To decouple the ni collinear sectors, we have to elimi- 
nate the products of usoft Wilson lines in cii^ . Using 

ni-nt _ . ni-fit , . 

ni - Aus = — — "-f^«s H 7. — "f^us + (ni)u-Aus 



2 2 
nt-Aus [l + 0{Xt)] , 



(3.39) 



10 



it follows that 



F„i (a;) = P exp 



15 



ds TifAusix + srit) 



and therefore 

YlY,,,=l + OiXt). 



0{\t) 
(3.40) 

(3.41) 



Using Eq. p.4ip in Eq. p.38p . the rii coUinear sector also 
decouples from the csoft one. We have now arrived at 
the same point as in Eq. p.24p . The remaining coupling 
of the (usoft-decoupled) csoft gluons to the ni coUinear 
sector via ni - An} is eliminated using the additional csoft 
field redefinition in Eq. p.25p . 

It is essential to perform the field redefinitions in this 
order. If we first perform a csoft field redefinition on ni, 
we would get a term 

Xl^gm-AusXn, (3.42) 

in Cm , which cannot be eliminated anymore by a usoft 
field redefinition. The fact that we have to perform the 
usoft field redefinition first and that it requires us to ex- 
pand in At, shows that this step is really linked to the 
SCET above the scale ^/t, which has At as its expansion 
parameter. 



IV. FACTORIZATION FOR e+e 



3 JETS 



In the previous section, we constructed a new effective 
theory, SCET+, which extends SCET with an additional 
mode that has csoft scaling. As a concrete example, in 
this section we apply SCET+ to 3-jet production in e+e" 
collisions. We are interested in the kinematic configura- 
tion shown in Fig. [5] We use A^-jettiness [6] with = 3 
to define the exclusive 3-jet final state, where the indi- 
vidual 3-jettiness contributions 71 of each jet determine 
the mass of the jets. We show how the factorization in 
SCET+ works and how the logarithms of the scales Q, t, 
and Ti are simultaneously resummed. 

We note that the applicability of SCET+ is not lim- 
ited to the class of A^-jettiness observables. However, 
A^-jettiness provides a convenient observable well suited 
for factorization because it is linear in momentum, does 
not depend on additional parameters (such as a jet ra- 
dius R), and covers all of phase space (i.e., there is no 
out-of-jet region). 



A. Definition of Observable and Power Counting 

1. Observables 

In terms of the lightlike jet reference momenta in 
Eq. (|2.ip . A^-jettiness is defined as [6| 



Ttv = y'niin{2(7, -pfe} , q: 



Qi 



(4.1) 




FIG. 2: The kinematic configuration of e"'"e~ — >■ 3 jets we 
consider. The boundaries of the jet regions determined by 
3-jettiness are illustrated by the dashed lines. Note that the 
location of the boundaries depends only on the jet reference 
momenta qi. 



where for convenience we defined the dimensionless ref- 
erence vectors qi. Tzv divides the phase space into A^ jet 
(or beam) regions, where a particle with momentum pk 
is in jet region i \i qi ■ pk < qj -pk for every j ^ i, i.e., the 
particle is closest to qi. The boundaries between the jet 
regions are illustrated by the dashed lines in Fig. [2j The 
Qi in Eq. (|4.ip are hard scales, such as the jet energies, 
Pt, or the total invariant mass. Different choices of Qi 
give different reference vectors qi, which lead to different 
choices of the distance measure used in dividing up the 
phase space into jet regions. The distance between two 
different jets is measured by the dimensionless quantity 



2q, ■ qj 



QiQj 



(4.2) 



For Qi ~ Ei, we have a geometric measure with qi = rii, 
and Sij — 2ni ■ nj measures the angle between jets i and 
J. For Qi — Q, we have an invariant-mass measure and 
Sij = Sij/Q"^ is equivalent to the invariant mass between 
jets i and j. By using qi and Sij we will keep our notation 
measure independent. (We will specify specific conditions 
on the used measure when necessary below.) 
We can write Tn as 



(4.3) 



where % is the contribution to Tzv from the ith jet region, 
which is given by 



%^2qi- '^Pk Qi{Pk) ■ 



(4.4) 



Here, the function 



Q^{p) = \{o{qo■p-q^■p) (4.5) 

imposes the phase space constraints for a particle with 
momentum p to lie in jet region i. Note that this con- 
straint only depends on the jet reference momenta (in 
addition to p itself) . 
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In the following, we will consider the cross section dif- 
ferential in each of the 71, as well as the minimum dijet 
invariant mass, t, and the jet energy fraction, z, defined 
as 



El 



El + E2 



(4.6) 



where the observed jets are numbered such that t = S12 
and El < i?2 ■ Since experimentally we cannot determine 
the type of hard parton initiating a jet, we will sum over 
all relevant partonic channels in the end. For simplicity, 
we also integrate over the three angles which together 
with t and z describe the full 3-body phase space of the 
three jets. (Two angles determine the overall orientation 
of the final state with respect to the beam axis. The 
third angle can be taken as the azimuthal angle of the 
two close jets.) 



which upon squaring yields 



(Qi + Q2)Q3 



(4.10) 



The dijet invariant masses S13 and S23 in the SCET above 
^/t are thus given by 



S13 = QiQ3Sq, S23=Q2Q3SQ 

which we can also write as 
si3 = xQ'^ , S23 = (1 - x)Q'^ , X 



Q1 + Q2 



(4.11) 



V- (4-12) 



In particular, for the geometric measure with Qi = Ei, 
we have x = z. Note that by counting all Qi ^ Q we in 
particular count Qi, 2/(93 ^ 1, which is necessary to have 
a consistent power expansion, such that 



2. Power Counting m SCET 

We consider the regime where all jets have similar en- 
ergies, such that Qi ^ Q and z ^ 1 — and take the dis- 
tance between jets 1 and 2 to be parametrically smaller 
than each of their distance to jet 3, such that 



■St 512 512 ^^2 
SQ Sl3 S23 * 

are all counted in the same way. 

3. Power Counting in SCET+ 



(4.13) 



St = S12 < S13 ^ S23 1 , 

i = S12 < si3 S23 ^ Q^ 



(4.7) 



corresponding to Eq. (|2.12p . 

To define the power expansion in our two-step match- 
ing procedure we now have to specify some power- 
counting properties of the distance measure. In the fol- 
lowing, we assume that we have chosen a measure such 
that the large components in qi and 172 are equal up to 
power corrections, such that S13 = S23 -I- 0(At). This is 
always the case for a geometric measure, where are 
effectively angles. For the invariant-mass measure, this 
is satisfied if the energies of jets 1 and 2 are equal up to 
power corrections. For measures where S13 and S23 differ 
by an amount of 0(1), the factorization still goes through 
but will have a somewhat different structure from what 
we will find below, and we leave the discussion of this 
case to future work. 

The power expansion of the SCET above the scale -\/t 
in terms of Af is defined by choosing a common reference 
vector qt for jets 1 and 2 in the direction of n^, such that 







^t 


= 9i[H 


-O(At)] =g2[l + 


0(At)], 




= nt 


• 93 


= Si3[l 


+ 0(At)] =S23[1 


+ 0(At)], 


Sit 


= 2(7t 


• 91 


= Si[H- 


'O(At)], 




S2t 


= 2gt 


• h 


= st[l4 


'0{\t)\. 


(4. 



The choice of qt is constrained by label momentum con- 
servation in SCET, 



(Q,o) = (gi + Q2)C + Q3g3^, 



(4.9) 



To setup the power expansion in SCET+ in A (or equiv- 
alently 77) , we first note that the invariant mass of the zth 
jet is given by 7] 



= P^^ QrTr [l + 0{%/Q^)\ 



(4.14) 



where Pi = "^j^Pk^iiPk) and so the invariant mass is 
determined by %. Hence, the condition Sij in 

Eq. (|2.12p . which requires the jet size to be small com- 
pared to the jet separation, corresponds to TiQi Sij. 
The power-counting parameters A^ = m^/Q^ and ij^ = 
m? /t are then determined by 



^-QX\ 



^-Qv' 

St 



-^^Qr]X. 

VSt 



(4.15) 



Note that to keep the power expansion in At ~ X/i] con- 
sistent, we still have to use the same vector qt (or n^) as 
in SCET to define the csoft modes in SCET+. This also 
applies when expanding the usoft measurement in At [see 
Eq. ((42T|) below]. 

All quantities related to the hard jet kinematics that 
enter in the final factorized cross section are uniquely 
determined in terms of the observables t and z by the 
label momentum conservation of the coUinear fields in 
SCET, , 



(Q,0) = Qi 9^ + ^2 92^ 



Q3 93 ■ 



(4.16) 



Recall that the large components of the coUinear fields 
in SCET_|_ are determined up to C'(A), which means we 
have to keep terms of 0(At) in Eq. (|4.16p . For example. 
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for the geometric measure where we choose Qi = Ei so 
qi — Ui, Eq. (|4.16p leads to 

0,^(1-^). (4.17) 

B. Factorization 

The SCET factorization theorem for the cross section 
fully differential in the 71 for equally separated jets was 
derived in Refs. [1, H]- The derivation in SCET+ follows 
the same logic, but we now have to take into account the 
presence of the new csoft modes. 

We first separate % for each i into its contributions 
from the collinear, csoft, and usoft sectors, 

T, = T^ + %'' + Tr, (4.18) 

where the individual contribution from different sectors 
are defined by restricting the sum over particles k in 
Eq. (|4.4p to a given sector. We will now determine the 
resulting measurement and phase space constraints for 
each sector. They are most easily obtained by expand- 
ing the full-theory measurement in Eq. (|4.4p using the 
appropriate momentum scaling of each mode. 

For a collinear mode in sector j with momentum p^, 
the distance to jet i, qi -pj, is by definition minimized for 
i = j, so 

0- (P ■) = % , (4.19) 

and therefore 

fcei-coU ^* 

where Pf is the total momentum in the i-coUinear sector, 
and (up to power corrections) si is the total invariant 
mass in the i-coUinear sector. Note that there are no 
phase space constraints from the jet boundaries in the 
collinear sectors, which leads to inclusive jet functions, 
J{si), in the factorization theorem. 

The division of the full measurement for the soft de- 
grees of freedom between the csoft and usoft sectors is 
more complicated and is illustrated in Fig. |3] in the 6-(j) 
plane. In Fig. EJa) we show the full measurement as de- 
termined by gi.2,3 and 0i,2.3- The three figures on the 
right show the various soft contributions which we will 
discuss next. 

To determine 7^"^, we write the 7V-jettiness measure 
for jet 1 for a usoft mode with momentum p„s in terms 
of the reference vectors q^ and qt, 

si3 . Sit ^ ^ ^ 

qi ■ Pus = — qt ■ Pus + — qs ■ Pus + [qi )Ut ■ Pus 
SQ SQ 

= qfPus[l + 0[Xt)\, (4.21) 



where we used the power counting in Eq. (|4.8p and the 
fact that all components of Pus have a common scaling. 
The same is true for 72 and q2, which means we can 
replace gi_2 by qt in the usoft contributions 7^"*. 

To determine the boundary between the jet regions 1 
and 2 we have to compare qi ■ Pus with 52 • Pus ■ For this 
comparison the subleading terms in Eq. ()4.2ip become 
relevant, so we have to be more careful. At the leading 
nontrivial order in the power counting this comparison 
only depends on the relative orientation of qi and (72 in 
the transverse direction. To have a simple way of writing 
the constraint, we can choose qt such that qi and (72 are 
back-to-back in the _Lt transverse plane and define the 
angle (ptipus) as the angle between qi and Pus in that 
plane. Then, pus is in region 1 for cos (f)t{pus) > and in 
region 2 for cos 4>t{pus) < 0. To summarize, we have 

/cG usoft 
A; £ usoft 

rr = 2q3- Y PkQTiPk), (4.22) 

/eg usoft 

where the boundaries are given by 

er{p) = o{q3-p-qfp)e[cosMp)], 
erip) = 0{q3-p~qfP)e[- cos Mp)] , 

er{p)=0{qfp-q3-p). (4.23) 

These are illustrated in Fig.lSjb). The hatching in regions 
1 and 2 denotes the fact that 77"! ^^'^ defined in terms of 
the common qt rather than their own qi or 92- 

Note that the standard 2-jet soft function depends on 
only two variables, whereas ours depends on three. How- 
ever, the only information about qi^2 that is retained in 
the usoft measurement is their collective direction, given 
by qt , and their relative orientation, given by (pt . In par- 
ticular, the usoft measurement contains no information 
about the angle between qi and (72, or equivalently st, 
at leading order in the power counting. This is in di- 
rect correspondence with the fact that the usoft modes 
only couple to the ni and n2 collinear sectors through 
a common Wilson line in the rit direction. Physically, 
the usoft modes are not energetic enough to resolve the 
difference between the rii and 712 directions. As a re- 
sult, the usoft function will only be sensitive to the scale 
Ti/ \/§Q ~ QX^ but not 7i/\/st ~ Qv^: which is consis- 
tent with our expectations from the physical picture as 
discussed in Sec. IIIBI 

The rif csoft modes are by definition collinear with jets 
1 and 2, so as with collinear modes their scaling implies 
that they are always closest to either qi or (j2. Hence, 
only the boundary between jets 1 and 2 remains, so the 
csoft phase space constraints are 

er(p) = ^(?2-P-9i-p), 
eTip) = oiqi-p-<i2-p), 

er(p)=0, (4.24) 
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(a) 71 
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(b) rr 
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(c) T-^' 



(d) 



FIG. 3: Graphical representation of the different measurement functions in the soft sectors in the 9-(l> plane for the geometric 
measure. The regions with different colors represent the phase space regions identified by the Oi(p), while the stars represent 
the directions of the dimensionless reference vectors qi used to calculate the observable. The full 3-jettiness measurement is 
shown on the left. The hatching on the right indicates a region where a different reference vector than on the left is used to 
compute the 3-jettiness observable. The contributions from the different hatched regions cancel on the right. 



and the csoft contributions to % are given by 



/cG csoft 



/c(^ csoft 



q-C 

I?. 



0. 



(4.25) 



The csoft measurement is illustrated in Fig. ^c). We 
now have only two different measurements, 7i and 72. 
In regions 1 and 2 they are computed with their proper 
reference vectors qi and (72 , reproducing the correct mea- 
surement in Fig.[3ja) for jets 1 and 2. At the same time, 
a different measurement is made in region 3, as indicated 
by the hatching. However, in region 3 the csoft modes 
are far away from nt, and so can only have usoft scaling 
there. Hence, the zero-bin subtraction of the csoft modes, 
which removes the double-counting with the usoft modes, 
will remove this region of phase space. 

Taking the usoft limit of Eq. using Eqs. (IT^ 

and (I4.23p . we obtain the csoft zero-bin contribution 



fcGcsoft— >-usoft 

rr° = 2g,. PkQT'iPk) 



fc£csoft—> usoft 



rr° = 0. 



(4.26) 



where the sum runs over all momenta in the csoft sector 
that actually have usoft scaling, and 



er°(p) = ^[cos^t(p)], 

0r°(p) -^?[- cosMp)] 



(4.27) 



The pictorial representation of this measurement is 
shown in Fig. ^d). As for the naive csoft, there are 
only two different measurements, but as indicated by the 
hatching in all regions the measurement is now performed 
with a different reference vector than the one used in the 
full 3-jettiness measurement. The complete csoft contri- 
bution is given by subtracting the zero-bin contributions 
in Eq. (IT^ from Eq. (g^l]). 



From Fig. |3] one can see how the total soft measure- 
ment in the full theory is reproduced by the combination 
of the usoft and csoft measurements. The zero-bin csoft 
measurement cancels both the csoft measurement in re- 
gion 3 made with a different reference vector than q^ and 
the usoft measurements in regions 1 and 2 made with dif- 
ferent reference vectors than qi and q2- The remaining 
csoft contribution in regions 1 and 2 and usoft contribu- 
tion in region 3 make up the correct measurement. To 
see this, consider the contribution of a generic soft gluon 
with momentum p to 7i. Summing up all its contribu- 
tions, we find 

(rr-rr° + rr)(p) 

^2qi -p 9{q2 -p - qi-p) -2qt ■ p 6[cos (j)t[p)\ 

+ 2qfp9{q3-p-qfP) d[cos <j)t (p)] 
= 2qi ■ plQiip) + 0{q2 -p-qi ■p)0{qi ■ P - qs ■ p)] 

- 2qt ■ pe[cos(l)t{p)]0{qt -p-qs-p) 
= Tiip)[l + OiXt)], (4.28) 

where Qi{p) is given in Eq. (j4.5l) . A similar equation is 
obtained for 72. For T3 we find 

{Tr-Tr' + TDip) 

= 2q3 -pQsip) + 2q3 ■ p[0{qt ■ P - q3 ■ p) - 03(p)] 
= T3{p)[l + 0{\t)\, (4.29) 

We will see this cancellation again explicitly in our one- 
loop calculation below. 

To formulate the measurement of 7i at the operator 
level, we define momentum operators which pick out the 
total momentum of all particles in each region according 
to Eqs. gll), (Oni) . (H^ . and 



Pi = y^Pfc8t(pfc) 



k 

E 

k 



PkoriPk) , 



E 



k 



PkQtiPk) , 

'^PkSriPk). (4.30) 



The differential cross section in 7i, 72, T3 in SCET_|_ 
is obtained from the forward scattering matrix element 
of the operator in Eq. p.35p . 



OpM3{Ti.T2,T3)O+\0). 



(4.31) 
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with the 3-jettiness measurement function 



(4.32) 



Using 71 = Tf^ + Tf" +Tj^'' from Eq. (|4T^ together with 
Eqs. (11201), (Hmi), (1123, and the momentum operators 
in Eq. (|4.30p . we can factorize the measurement function, 



Mi{TuT2,%) 



4=1,2 



S3 



'3 ' 



(4.33) 



where the cohinear, csoft, and usoft measurement functions are 

1=1,2 

-T^^.^ ^ _ • Pr) n '^('^"' - ■ ■ (4-34) 

1=1,2 

This factorization of the measurement function together with the factorization of the operator discussed in 
Sec. nil CI allows us to factorize Eq. (|4.3ip into separate collinear, csoft, and usoft matrix elements. This is the 
cornerstone in obtaining the factorization theorem for the differential cross section. The derivation of the final 
factorization formula now only requires one to properly deal with the phase space sums over label and residual 
momentum and to provide an operator definition of all components in the factorization theorem. The required steps 
in SCET+ are straightforward and the same as in SCET, see Refs. [1, [IB [3, H^l ■ The final factorized cross section, 
differential in the %, t, and z is given by 



da 



dridr2dr3dtdz 



= ^E^2(Q',A*)^^(i,^,M)n / ds.J«,(s„M) 

X J dhdk2S'l{h,k2,fl)S2[Tl-^^-ki,T2 



S2 , ^ S-i 
W2 W3 



(4.35) 



Here, co — {4:Tralj^/3Q^)Nc J2q Qq ^^'^ tree-level cross section for e+e^ — ,> hadrons. 

Since jets initiated by different types of partons are not distinguished experimentally, we sum over the relevant 
partonic channels to produce the observed jets, which are labeled such that the minimum dijet invariant mass t is 
Si2 and El < £'2- The sum over partonic channels is denoted by the sum over k = {ki, K2, K3}, which runs over the 
four partonic channels k = {q, g, q}, {g, q, q}, {q, g, q} and {g, q, q}. For the first two channels, jets 1 and 2 effectively 
arise from a. q^ qg splitting, and for the last two from a. q ^ qg splitting. For each splitting there are two channels, 
depending on whether the gluon or (anti)quark has the larger energy fraction. (The contribution where the quark and 
antiquark form the two jets with the smallest invariant mass does not enter in the sum because it is power suppressed.) 

The hard function H2 is the squared Wilson coefficient of O2 from matching QCD onto SCET, and in our case is 
independent of k. The hard function H'^ is the squared Wilson coefficient of from matching SCET onto SCET+. 
The Jr. (s,^) are the standard inclusive jet functions in SCET and the soft functions S2 and denote the matrix 
elements of the usoft and csoft fields, respectively. 



(rr,rr,M) 



1 



NcCf 



(4.36) 



The soft functions implicitly depend on the reference vec- 
tors qi through the combinations sq and st, respectively, 
which is suppressed in our notation. The definition for S2 
is given for rit and 713 corresponding to a quark and anti- 
quark, respectively, but 5*2 itself is independent of k, i.e.. 



it is the same for (? O g, which only switches F o y^. 
The definition of 5+ is given for k = {g, g} for which 
Til = Uq, n2 = g and rit corresponds to a quark. The 
definitions for the other channels follow from the obvious 
interchanges of the appropriate Wilson lines. 
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In the next section we discuss all the ingredients in 
Eq. (|4.35p in detail, and obtain their explicit one-loop 
expressions. We also discuss the relation of the hard and 
soft functions in Eq. (|4.35p to the 3-jet hard and soft func- 
tions in SCET, and derive the structure of the anomalous 
dimensions of iJ+ and 5+ to all orders in perturbation 
theory. Readers not interested in these details can skip 
to Sec. IVII where we give the generalization of Eq. (j4.35p 
to pp ^ N jets or to Sec. IVIII where we present explicit 
numerical results for the dijet invariant-mass spectrum 
resulting from Eq. (|i35|) at NLL'. 



V. PERTURBATIVE RESULTS FOR e+e" 3 
JETS 

To exhibit the color structure and be able to easily 
generalize our results in Sec. IVII we will use the standard 
color-charge notation, where denotes the color charge 
of the ith external parton when coupling to a gluon with 
color A. In general • = T^T^ are matrices in 
the color space of the external partons. In particular 



Tf = 1 C, where 



Cq — Cq — Cp . 



Cg = CA- (5.1) 



In the following, we will have three external partons, qqg, 
for which the color space is still one-dimensional and the 
color matrices reduce to numbers. 



1, 
Ca 
2 



rp2 



T T 



^T^2 

g ~ 
2 



Ca ■ 



(5.2) 



Hard EYinctions 



As discussed in Sees. Ill 61 and IhTC] for the jet config- 
uration in Fig. [5] we are interested in, the matching onto 
the operator proceeds in two steps. This allows the 
dependence on the two parametrically different scales, 
t = Si2 <^ Si3 ^ S23 ^ to be separated. 

In the first step we match at the hard scale Q from 
QCD onto SCET as shown in Fig. 1(b) In our case we 
match the QCD current, ■07'^^/', onto the SCET two-jet 
operator, O2 = Xntl'^Xnq, by computing and comparing 
the qq matrix elements in both theories, 

Xqcd(0^(/^) =C2(Q^M)(99|O2(M)|0). (5.3) 

Here, 02{n) denotes the MS renormahzed operator. This 
matching is well-known (see, e.g., Ref. Q for a de- 
tailed discussion) and was first performed at one loop 
in Refs. [l^ll^. The resulting matching coefficient is 



and satisfies the RGE 

/^^c•2(Q^M)- 

The one-loop anomalous dimension is given by 



7c..(Q^M)C2(g^A^). 



7C2(Q^A^) 



,{^^)c, 

47r 



41n- 



-g' - io 



6 



(5.5) 



(5.6) 



The hard function H2 [Q^ , jj) in Eq. (|4.35p and its anoma- 
lous dimension are given by 



H2{Q^^i) = \C2{Q^^i)\'' , 

7ff,(Q",M) =2Rc[7c2(Q',A*)] 



(5.7) 



We then run down to the scale ^/t, and match from 
O2 in SCET to the operator in SCET+, as shown 



in Fig. 1(b) In principle, this matching is computed in 



an analogous way by calculating the relevant 3-parton 
matrix elements in both theories (suppressing any spin 
indices), 

{qq9\02{^im = C'l{t,x,fi) {qqg\O+{fi)\0) . (5.8) 

The full one-loop calculation for the matrix element of O2 
is quite involved. However, we can extract the one-loop 
result for the hard function H'^, given by 



(5.9) 



using the known one-loop result for e+e^ — > 3 jets from 
Ref. [131 • Since the operator matching Eq. (|5.3p is inde- 
pendent of the final state, it follows that'^ 



\Mqcb{0 qqg)^ 

^H2{Q^^i)\{qqg\02{^im{' 
^H2iQ'\fi)H^it,x,fi)\{qqg\O+{ii)\0)\\ 



S12<SS13~S23 

-l2 ,,\ |/„-„|/n f ,.\\n\\2 



(5.10) 



In pure dimensional regularization, the virtual one-loop 
corrections to the bare matrix element are scaleless 
and vanish. Hence, the renormahzed matrix element of 
0^(/i) on the right-hand side is given by the tree-level 
result plus the counter-term contribution, which effec- 
tively supplies the proper 1/e divergences to cancel the 
IR divergences in the left-hand side matrix element. The 
remaining finite terms then determine the one-loop cor- 
rections to H'^{^). 

For K = {q,g,q}, we take t = Sqg, S13 = s,,, 
•S23 = Sqg. Expanding the one- loop virtual corrections 
for |A^qcd(0 ^- qqg)\'^ from Ref. 27] in the limit t ^ Q'^ 
with Sqq = xQ^ and Sqg = (1 — x)Q'^ [see Eq. (|4.12|) ]. 
and combining them with the one-loop corrections to 
H2{Q'^,fi), we find 



C2{Q^^i) = l + 



as{^i) Cf 



4tt 



3 In 



In^ 



(^) 



(^) 



-8 



21 



(5.4) 



^ Since we are only interested in the cross section integrated over 
angles, we can consider the spin-averaged matrix element, which 
removes the dependence on the azimuthal angle in the q ^ qg 
splitting. Including this dependence requires to explicitly take 
into account the spin structure of Oj^. 
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H 



{9,3,9} 



2 



2tt 



t l-x 



27r 



V 2 



Ci.) (^2 In In x + In^ x + 2Li2(l - x) 



Ca 



2 V 



21n- 



■ ln(l -x)+ ln^(l - x) + 2Li2(x) + (Ca - Cj.) 



l-x 

1 + a;2 



(5.11) 



The overall factor of here is included to make iJ+ di- 
mensionless. Note that at tree level takes the form of 
the common q qg splitting function. As discussed be- 
low Eq. (j6.18p . beyond tree level it is related to universal 
splitting amplitudes (which are not the same as splitting 
functions). The results for the other partonic channels 
are given by 



H 



+ 

{9,9,9} 



(5.12) 



{t, X, = hI^'^'"^ it, X, m) = i^l"'^'"^ (t, l-x,y). 



{<i,g,a} I 



They follow from the fact that |A^qcd(0 qqg)]^ is 
symmetric under the interchange q q. 
The hard function satisfies the RGE 



fi—H^{t,x, ^l) = -f^^{t,x, ^l) H^it^x, ^l) . (5.13) 

At one loop we find for k — {q, g, q} using Eq. (|5.1ip 



2tt 



t 



2Ca In — + 4 (Cf - In x 



Ca\ 



+ 2CaH1 - x) ~ l3o 



(5.14) 



where f3o = (IICa— 4rFn/)/3. For general k, the anoma- 
lous dimension can be written as 



lH+ (t, X, m) 



2tt 



4Ti • T2 ln4r 4-4Ti-T3 Inx 



4T2-T3 Hl-x)+po 



(5.15) 



Its all-order structure is derived from consistency in 
Sec.lVDl 



B. Jet Functions 

The jet functions are given by the matrix elements of 
coUinear fields, and are the standard inclusive jet func- 
tions as in many other SCET applications. We give the 
one-loop renormalized jet function in MS for complete- 
ness 



An 



3 p 



Jg{s,^l) = S{S) + 



47r 



ml 



(5.16) 



where /3o = (IIC^ — 4:Tpnf)/3 and Cn{x) denotes the 
standard plus distribution, 



£„(x) = 



(x) ln"x 



The jet functions satisfy the RGE 

with the anomalous dimensions 

Yj{s,f^) = -2Cjrcusp[as(M)] -J^o(^) 



(5.17) 



(5.18) 



(5.19) 



where rcusp[Q!s] is the universal cusp anomalous dimen- 



321 1 given in Eq. (jCSp . and the noncusp terms are 



given in Eq. (|C13|) . 



Soft Functions 



The usoft and csoft functions describe the contribu- 
tions to the observable from particles softer than the jet 
energies. Unlike collinear modes which contribute to only 
a single jet, the soft modes can contribute to all jets. This 
means that these modes are sensitive to the invariant 
masses between jets. The csoft modes, while having 
smaller energy than the collinear modes, have collinear 
scaling and are needed to describe the soft interactions 
between the nearby pair of jets, because the usoft modes 
have too small energy to resolve these two jets. The 
results for the soft functions can be written for general 
AC, i.e., without having to specify a particular channel, 
since the dependence on the parton species solely arises 
through the SU(3) color representations of the Wilson 
lines. 



1. The Ultra-Soft Function S2 

The operator definition of the usoft function ^2 is given 
in Eq. (I4.36|) . with the measurement function given in 
Eq. (j4.34p . At one loop, the relevant integral we have to 
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ni 



®: 




■n2 



n2 



ni 





(a) 



(b) 



(c) 



nt 




nt 



rii 



:® 



(d) 



FIG. 4: One-loop diagrams for 5'+. The vertical line denotes the final-state cut. There are also diagrams analogous to (c) and 
(d) with the gluon coupling to the n2 Wilson lines. Virtual diagrams and diagrams with the gluon coupling to the same Wilson 
line are not shown. 



compute is 



2(T,-,T.).T3(^)V/^2.5bW) 



"-3 ■ nt 



(na ■p){nt ■ p) 



Mr{ki,k2,k3) 



(5.20) 



For the color factor we have used that with respect to 
the external 3-parton color space, the total color charge 
carried by the nt Wilson line is Ti-f-T2, i.e. the combined 
color of partons 1 and 2. 

The usoft region for jets 1 and 2 is determined by 
Eq. (|4.23p . where the boundary between jet 3 and jets 
1 and 2 depends only on qt and (73 . The division of the 
combined nt region between jets 1 and 2 is given by the 



additional 9[± cos (f>t{p)], whose only effect is to divide 
the azimuthal integral in half. In d = 4 — 2e dimensions 
one gets 



■7T/2 



ht sin 0t 



h sin 2' <l)t 



7r/2 

1 r 
2./0 



/ d0t sin"2^ (f>t . (5.21) 
Jo 



Hence, the nt hemisphere contribution is split in half 
between jets 1 and 2. 

The final result for the renormalized usoft function at 
NLO is 



S2{ki,k2,k3, /x) = 6{ki) S{k2) Sih) + (Ti + T2) • T3 



An 



ki 



Sik2)S{k3) 



<5(fci) ,5(fc3) + £1 (^^)(5(fci) 5{k2) 



-y<5(fcl)<5(fc2)5(fc3)J>, 



(5.22) 



which is simply the sum of two hemisphere contributions. We can see explicitly that S2 depends only on the scale 
ki I V^sq ■ It satisfies the RGE 



^d/I "^^e^i'^^a, k'i.ii) = y d/c'i dfcadfcg 752(^1 - k[,k2 - k'2,k3 - fcg. 



^j.)S2{k[,k'2,k'3,fi) , 



where the anomalous dimension at one loop is given by 

75.(fcl,fc2,fc3,A*) = -^(Ti+T2)-T3 ^ 



2Co(^^)siki)6{k2) 



lco(^^)5ik2)5(k,) + Co(^^)s{k,)S{k,) 



(5.23) 



(5.24) 



2. The Collinear-Soft Function 



The definition of the csoft function 5+ in terms of a matrix element of Wilson lines is given in Eq. (j4.36p , with the 
measurement function given in Eq. (I4.34[) . The calculation for 5+ is more nontrivial due to additional csoft Wilson 
lines V„^ , and we therefore provide some more details. 
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There are two basic types of diagrams at one loop, shown in Fig. |4l In the diagrams shown in Figs. 4(a) and 
|4(b)[ a gluon is exchanged between the X Wilson lines in the rii and 712 directions, which corresponds to a csoft 
gluon exchanged between the nearby jets. These diagrams are the same as in a usual soft-function calculation. The 
analogous virtual diagrams vanish in pure dimensional regularization, and the diagrams with the gluon attaching to 
the same Wilson line vanish due to — 0. The diagrams do not require a zero-bin subtraction, and their contribution 
to the one-loop renormalized csoft function is the usual hemisphere contribution 



47r 



Ti T2 



-A I 



Similarly, the contribution to the anomalous dimension from this diagram is the hemisphere contribution. 



7g+42(fci.fc2,/i)--^Ti.T2 



Co 



k_2 



(5.25) 



(5.26) 



The second type of diagram comes from exchange of a gluon between the V and X Wilson lines, as shown in Figs. 4(c) 
and 4(d)[ These diagrams have a nontrivial usoft limit, which means we must perform a zero-bin subtraction to remove 
double counting. As discussed in Sec. IIIH the zero-bin limit is obtained by expanding ni in terms of nt, with the 
zero-bin measurement obtained from Eq. (j4.26p . We focus on gluon exchange between Vr^ and X„j ; the results for 
changing ni — )■ n2 are analogous. Subtracting the zero-bin contribution from the naive part of the diagram yields 



-2 Ti • T, 



in 



nt ■ ni 



{nt ■p){ni -p) 



W{k^M) 



rit ■ nt 



{nt ■p){nt -p) 



(5.27) 



Here, Ti is the color charge carried by X^- Since S+ is diagonal in color, the color charge carried by Vnt is 
— Ti — T2 = T3. The result for Eq. (|5.27p can be extracted using the results of Rcf. [7] in the limit St sq. We split 
up the phase space for the naive and zero-bin csoft contribution into regions qs ■ p > qi ■ p and ■ p < qi ■ p, where in 
the latter region the naive and zero-bin contributions cancel. In the region qs ■ p > qi ■ p the naive csoft contribution 
is given by the sum of the hemisphere and nonhemisphere contributions iSjg'j^pjjj; -I- >S'{3''2 of Ref. Q expanded in the 
limit St <C SQ. It is straightforward to calculate the zero-bin contribution in Eq. (|5.27|) for q^ ■ p > qi ■ p. Taking the 
difference between these terms gives the total contribution to the renormalized one-loop function from V^n and 
Xn-i exchange 



S'+llk^i'^a,^) = 



Qs(m) 
47r 



Ti -Ta 



^/:i(^)<5(fc2) - -^cA-^)5{k,) - '^Sik,)S{k2) 



27^ 

3 



The corresponding contribution to the anomalous dimension is 



ki 



Co{j]Sik2)-Co{^]5{ki 



(5.28) 



(5.29) 



where ^ is a dimension-one dummy variable which is only needed to make the argument of Co dimensionless, but 
cancels between the two terms. The analogous contribution with a gluon exchanged between Vm and the 712 Wilson 
line is the same with the replacement 1 o 2. 

Combining everything, the final result for the one-loop renormalized function becomes 



Slik,,k2,^^)^^ik,)s{k2) + 

+ (Ti - T2) 



«^(m) 
An 

•T3 



Ti T2 



-A I 



ki 



st^^ 



ki 



5{k2)~Ci(^ 



)5{k2) + 
k2 



-A 



st^J. 

5{ki) 



^\/ StU'' 



k 

(Ti 



n' 
' "3 
27r2 
~3~ 



5{ki)5{k2) 

S{ki)S{k2) 



The RGE for 5*+ has the form 



We can see explicitly that 5*+ depends only on the scales ki / Vst ■ 

^dfei dk'2 7L {kl-klk2~k'2,fi) SI {k[ ,k'2,fi), 



M^5'+(A:i,fc2,M) 
where the one-loop anomalous dimension is given by 



7S (fci,fc2,A*) = -- 



in 



Ti •T2 



Cr 



ki 



6{k2 



Cr 



S{ki) 



(5.30) 



(5.31) 



(5.32) 



(Ti - T2) • T3 - 



Co{j)sik2) 



C 
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3. Soft Functions with Single Argument 



Hard-Function Consistency and Derivation of 70^ 



For our numerical analysis in Sec. I VIII we project the 
soft functions onto the sum of their arguments, 



S'2(fc,A*) = J dkidk2dk3 S2{ki,k2,k3, ^) 

X (5(fc- fci - fca - fcs), (5.33) 
S+{k,^) ^ / dfcidfc2 S1{ki, fc2,^) S{k - fci - ^2) . 



From Eqs. (lO^ and (jOg]) . we obtain their NLO ex- 
pressions, 



S2ik, ^i) ^ s{k) + 



as{^^)CF 
An 

16 . / fc 



SQfiy 3 



S+ik,,,)=S{k) + ^4^{~CA 



An 



-A 



^-^Sik) 



+ Cp^S{k) 



(5.34) 



Note that this projection removes the dependence on 
(Ti — T2) • T3, which makes S+{k,fi) independent of 
K. The single- argument soft functions satisfy the RGE 



M^'S'(fc, /^) = y dfc' 7s(/c - k\ n) S{k', n) , (5.35) 

where the anomalous dimensions after projecting onto k 
simplify to 

Q!s(m)Cf 16 / A: \ 
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7s+(fc,Ai) 



as(/i)CA 8 



An \/TtiJ, V 



(5.36) 



D. All-Order Anomalous Dimensions 

In this section we discuss the consistency constraints 
on our factorized cross section in Eq. (|4.35p . This allows 
us to derive the general form of the anomalous dimen- 
sions for the SCET+ matching coefficient, C+, and csoft 
function, 5+ , which are the new ingredients in the factor- 
ization from SCET+ . In particular, we demonstrate that 
the convolution of the csoft and usoft functions at one 
loop reproduces the known result for the 3-jettiness soft 
function in regular SCET in the limit S12 ^ S13 ^ S23. 
This demonstrates that the csoft modes are necessary for 
SCET+ to reproduce the correct IR structure of QCD in 
this limit. We then show that the factorized cross section 
obeys exact renormalization group consistency. 



The factorized 3-jettiness cross section in SCET is 
given by [71] 

d(T 



dridr2d7i dtdz 



n2 



Here, all dijet invariant masses are counted as ~ Q^. 
This means that the hard function, H^, is evaluated at 
their exact values given in terms of t and z, 

si2 = t, si3 = zQ^ - (1 - z) t , S23 = (1 ~ z)(3^ - 2i , 

(5.38) 

which follow from momentum conservation for e~^e~ 3 
massless jets. At tree level. 



H. 



(Sl2, Sl3, S23,At) 

asipjCp (si3 -I- 523)^ + (si2 + 513)^ 



2n 



S12S23 



(5.39) 



In SCET, all loop diagrams contributing to the bare 
matrix element of {qgq\O3\0) vanish in pure dimensional 
regularization, and consequently the 3-jet hard function 
in SCET, i?3({sij}, ^), is directly given by the IR finite 
terms of the fuU QCD amplitude |A^qcd(0 -> qqg)\^- 
Comparing with Eq. (I5.10p . it follows that the hard func- 
tions in SCET and SCET+ to all orders in perturbation 
theory have to satisfy 



(5.40) 

At tree level, this can be seen immediately: to expand 
Eq. (I5.39|) in the limit S12 ^ S13 ^ S23, we set S13 = xQ"^ , 
523 = (1 — x)Q'^ [see Eq. (I4.12p ]. and t = S12 and drop 
any terms subleading in t/Q^, which gives the tree-level 

resuh for i/j,''^'^^ (i, x, ^i) in Eq. (1011) . 

The above argument also applies directly to the Wilson 
coefficients before squaring them, so 



Sl2^Si3~S23 



C2{Q^^Ji)Cl{t,x,^Ji). 



(5.41) 

Taking the derivative with respect to /i, it follows that 



Sl2^Si3~S23 



7C2(Q^A^) +7£, (i,^;,M) 



(5.42) 

The general all-order forms of the anomalous dimensions 
7C2 and 7^3 are [H,!!!,!!! 

7C2(Q',m) - -rcuspK(M)] (T1+T2) -Ta ln ~^'~'° 



+ 7cK(Ai)] +7cK(m)] 
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-i<J 

+ 7c ["^ (m)] + 7c ["s (a*)] + 7c ["^ (a*)] > 

(5.43) 

where the individual quark and gluon contributions in 
the noncusp terms are given in Eq. (jC14p . Compared to 
Eq. (|5.6p we have identified the color structure in 7C2 as 



T, • T, 



gg 



(Ti+T2)-Ta 



ggg 



(5.44) 



Here denotes the combined color charge of the quark 
or antiquark that splits into partons 1 and 2, and TfT^ is 
evaluated in the corresponding 2-parton qq color space, 
i.e., Tq ■ Tqlqq = —Cp- lu thc sccoud step, we wrote 
the same total color charge using the individual color 
charges of the daughter partons 1 and 2, which are now 
evaluated with respect to the 3-parton qqg color space. 
Exphcitly, using Eq. (|5.2p we have (T, + Tg) ■Tq\qqg — 
[Ca/'^ — Cp) — CaI"^ = —Cp, and with the same result 
for q q. 

Using Eqs. (|5.42p and (|5.43p and expanding 7^,^, we 
obtain the general form of valid to all orders in 

perturbation theory. 



7C+ 2;, /i) 



[a,(Ai)]Ti •T2ln■ 



-^- iO 
1^ 



+ 7c+K(m)'2;] , 
7£ [a„ x] = -reusp[a.] [Ti -Tg Inx + T2 -Tg ln(l - x) 



+ 7&[«s] 



(5.45) 



Note that this provides a nontrivial example of a hard 
anomalous dimension, where the nonlogarithmic term, 
7c_|_ [oisj x\, depends on a kinematic variable, whose over- 
all coefficient however is still determined by Fcusp- At 
one loop, Eq. (I5.45[) reproduces Eq. (|5.15p exactly using 
that 7^_^(t,a:,/i) = 2Re[7^^ (i, x, /z)]. 



2. Soft-Function Consistency and Derivation of 75^ 

In Sees, [n] and IIIII we have seen that SCET+ arises 
from expanding SCET in the limit t ^ Q. It follows 
that the SCET+ 3-jet cross section in Eq. (I4.35P has to 
reproduce the 3-jet cross section Eq. (|5.37p computed in 



SCET when the latter is expanded in the limit S12 <ti 

Sl3 S23, 



da 



SCET 



= dcr 



SCET4 



S12>SS13~S23 



(5.46) 



(This is exactly analogous to the statement that the 
SCET cross section must reproduce the QCD cross sec- 
tion expanded in the limit m <^ Q.) As we have seen 
above, the product of hard functions in SCET_|_ repro- 
duces the full SCET hard function, and the jet functions 
are the same in both cases. Hence, for the cross sections 
to satisfy Eq. (I5.46p . the soft functions have to satisfy 



SfCsl3=S23 

dk[dk'2S2{ki-k[,k2 



(5.47) 

k',,k3,fl)S^{k[,k'^,fl). 



For the soft functions the limit S12 <^ S13 ~ S23 is taken 
using Eq. (|4.8p by setting st = S12, S13 = S23 = sq and 
expanding in Sf <C sg. 

The fact that the hard and soft functions separately 
factorize in the limit t ^ Q as in Eqs. (15.411) and (I5.47P 
is a direct consequence of factorization in SCET and 
SCET+. Since the soft sectors in both theories are decou- 
pled from the coUinear sectors, the soft sector of SCET+ 
has to reproduce the soft sector of SCET expanded in 
i <C Q. Since the factorization applies also in the kine- 
matic region where the soft functions become nonpertur- 
bative, the relation in Eq. (|5.47p between the soft func- 
tions in the two theories holds both at the perturbative 
and also the nonperturbative level. 

We can check explicitly that Eq. (|5.47l) is satisfied by 
our one-loop results. Since SCET correctly reproduces 
the IR structure of QCD, this also provides an explicit 
demonstration at the one-loop level that the csoft modes 
are necessary to reproduce the IR structure of QCD in 
the limit m <^ t <^ Q, and thus that SCET+ is the 
appropriate effective theory of QCD in this limit. 

The full A^-jettiness soft function at NLO has been 
calculated explicitly in Ref. where the final result is 
given in terms of a single integral, which can be evalu- 
ated numerically. In Ref. [s^ a general algorithm was 
developed to calculate a wide class of soft functions with 
an arbitrary number of collinear directions numerically. 
In the limit st <C sq, the required integrals for in 
Ref. 0] can be obtained analytically, and we find 



S'a (fci, fc2,3 ,/^) 



St^Sl3=S23 



5(fcl)J(fc2)<5(fc3) 



Ti • T3 



An 



fcl 



Ti-T2 

6{k^) 



ki 



V(fc2) + -^cJ^)6{k,) - \5{k,)5{k2) 

S{k2) 



Sik,) 



S{k,)-—6{k,)Sik3) 
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cJ^-)sih) + -^/:i(-^)<5(fc2) - ^(5(fc2)<5(fc3) 



(Ti ^Ta) •T3 21n(^)i 



/:o(y)<5(fc2)-/:o(^)^(fci) 



(Ti + T2) • T3 5(fc2) 5(fc3) !> . (5.4 



Here, the first three terms proportional to • Tj are the hemisphere contributions, which contain the explicit /i 
dependence. The last two terms come from the nonhemisphere contributions, where the dummy variable ^ again 
cancels between the two terms and is only needed to make the argument of £n dimensionless. As a cross check, 
we have compared this result with the numerical result obtained using Ref. [35|, and the two agree in the limit 
St ^ S13 = S23 = SQ. As expected, this soft function depends on both st and sg, and there is no choice of 
renormalization scale for which the logarithms of St/sq are absent. 

Since the soft functions at tree level are all just 5 functions, Eq. (|5.47p simplifies at one loop to 



5*2 (fcl, ^2, k3,^l) + 57 '{ki, fc2,Ai)(5(fc3) . 



(5.49) 



Subtracting the as corrections in Eqs. (15.301) and (I5.22[) from those in Eq. (|5.48p . the terms proportional to Ti • T2 
and those involving £1(^3) or only S functions immediately cancel. For the remaining terms, we obtain 
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(Ti -T2)-T. 



4 



f^)sik,)-c,(^)s{k,) 



k-2 



fcl 

's't^J, 



ko 



21n(|)i[£o(f)^(M-^o(f)^(fci) 



cA-^)6{k,)-cA^)6{k,) 



= 0. 



(5.50) 



where we used the rescaling identity 



XCi{\x) = Ci{x) + \n \ Ca{x) + -lnXd{x) 



(5.51) 



Thus, Eq. (|5.47p is satisfied by our one-loop results. 

We can also use Eq. (|5.47p to derive the all-order structure of the anomalous dimension of S+ . Taking the derivative 
of Eq. ()5.47p with respect to fi, we get 



753(^1-^2,^3,^^) 



St<Ksi3=S23 



752(^1, ^2,^3,^) +7s+(^l'^2,M)<5(fc3) 



(5.52) 



The all-order structure of the anomalous dimension of the A^-jettiness soft function was derived in Ref. Q- For 
in this limit we have 



7S3(fcl,*2,fc3,M) 



St<Ssl3=S23 



-2rcuspK(M)]<^Ti •T2 



1 



+ Ti • T3 



T, • T3 



1 



Co 



k2 



J^)6ik2)+Co(^)6{k,) 



Co(^-)s{k,) + Co(^^)6{k,) 



Siks) + Co 



-)S{k2 



ka 

5(k2) 

5(ki) 



Sik,) 



+ -fl[aM]Kki)Sik2)Sik3). 



(5.53) 



To deduce the all-order structure for 752, we first note that upon projecting S2{ki, ^2, ^3) onto kt = ki + ^2, 



S2{kt,k3,fj,) = J dki dk2 52(^1,^2, k3,^)S{kt - fci - fc2) , 
it reduces to the normal 2-jettiness soft function. Therefore, we know that to all orders 



(5.54) 



dfci dfc2 752(^1, fc2, ^3, fi) 6{kt - ki - k2) 

1 



= -2reusp[as(M)] (Ti+T2)-T3 



Cn 



J^^)s{k,) + Co(^^)sih) 



--1sAasm5{kt)Kk^)- (5.55) 
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From its definition, we know that the full 6*2(^1, ^2, ^3) is symmetric in ki and k2, and since the distinction between 
ki and k2 only comes from the measurement function, this symmetry cannot be changed by the renormalization and 
the anomalous dimension must therefore also be symmetric in ki and fc2 • Furthermore, from Eq. ()5.52p we know that 
the dependence on ^3 must exactly cancel between 75^ and 75^. The most general form of 75^ that satisfies these 
requirements, Eq. (|5.55p . and is only single-logarithmic in /i is 



75.(fcl,fc2,fc3,M) --2rcuspK(M)] (Ti+T2)-T3^^|i£o(^^)'5(fc2)'J(fc3 



Co 



S{ki) S{k2) \ + 7s. S{ki) S{k2) Siks) 



lco(^-)s{k,)Sik,) 

(5.56) 



The noncusp terms in Eqs. (|5.53p and (|5.56l) are 

^U^s]^- ^ (7}K]+27^K]) =0 + O(a2), 75,[a,] ^ (7}[a,]+27^[a,]) =0 + O(a2), (5.57) 

where 7j[as] and 7p[as] are the noncusp terms in the jet and hard anomalous dimensions Eqs. (|5.19p and (I5.43p . and 
are given in Eqs. (IC13I) and (IC14I) . This form of Eq. (|5.56l) agrees with our one- loop result in Eq. (j5.24p . Taking the 
difference between Eqs. (I5.53P and (|5.56p . we obtain the general form of 7^^, 



lsAh,k2,^l) = -2reusp[a.(Ai)]Ti -Ta 



1 



7S [as , ki , k2] = -Feusp [as] (Ti - T2) • T 



1 

1 



Co 



5{k2) + Co 



5{ki 



+ 7s+K(M),fci,fc2], (5.58) 



Co{j)5{k2)-Co{j)5{k, 



-(75M + 27S[a.])<5(fci)<5(fc2), 



which again agrees with our explicit one-loop result in 
Eq. (|5.32l) . The part of the anomalous dimension of 
5"+ which does not explicitly depend on /i has a more 
complicated structure than for 5*2 and S'3. It has a 
nontrivial color structure and dependence on the kine- 
matic variables fci and fc2j which effectively behaves as 
ln(fci/fc2). The coefficient of that dependence is how- 
ever still determined by Fcusp- This is the soft analog 
of what we saw for the anomalous dimension of C+ in 
Eq. (|5.45p . which contains terms like \^{si^/Q^) = ln(a;) 
and ln(s23/(5^) ln(l - x). 



3. Combined Consistency of Factorized Cross Section 



which is derived by taking the derivative of Eq. (I4.35P 
with respect to ^, and following the same steps as in 
Ref. to derive the analogous relation for the cross sec- 
tion in SCET given in Eq. (|07)) . 

Since the SCET cross section satisfies the RGE consis- 
tency, we already know that Eq. (|5.59p must be satisfied 
as well. Nevertheless it is an instructive and straightfor- 
ward exercise to show that Eq. (|5.59p is indeed satisfied 
to all orders by the results for the anomalous dimensions 
given in Eqs. (EH]), (|535)) . and (1^351) . 

The cancellation of the different logarithmic dependence 
in the hard, jet, and soft functions for the three color 
structures Ti • T2, Ti • T3, and T2 • T3, now happens as 
follows, 



As we have seen above, the sum of the hard and soft 
anomalous dimensions in SCET+ each reproduce the 
hard and soft anomalous dimension "fc2 ^-nd 7S3 in SCET 
in the limit S12 ^ S13 ~ S23- The full cross section in 
Eq. (|4.35p is a physical observable and cannot depend on 
the arbitrary renormalization scale /i. This implies that 
the anomalous dimensions must satisfy the consistency 
relation 

= 2Re [7c, (Q^ m) + lc+ (t, A*)] H '^(^') 

i 

+ IS2 {ki,k2,k3, n) + 75^ (fci, fc2, m) S{k3), (5.59) 
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The \n{Qi£^/fi^) terms are supplied by the jet functions. 
Note that this cancellation crucially relies on a consistent 
power expansion in Xt, as in Eqs. (|4.8p and (|4.1ip . which 
implies S13 = xQ" ^ = Q 1Q3SQ and S23 = (1 - x)Q^ = 
Q2Q3SQ, so Eq. (|5.60p is satisfied exactly without re- 
quiring any further expansion. 
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VI. GENERALIZATION TO pp ^ iV JETS 

In Sees. IIVI and IVl we have applied our new effective 
theory SCET_|. to the simple case of e+e" — > 3 jets. This 
allowed us to discuss in detail how SCET_|_ is applied to 
derive the factorized cross section, and to obtain all its 
ingredients at NLO. In this section, we extend our dis- 
cussion to the general case of pp — )■ TV jets plus additional 
leptons relevant for the LHC. In particular this requires 
adding hadrons to the initial state, as well as generalizing 
to more final state jets and the resulting more compli- 
cated color structure. 

The key ingredients needed to derive the factorization 
theorem are the same here as in our 3-jet example in 
Sec. HVl We have to define a consistent power counting, 
determine the relevant operators, and show the factoriza- 
tion of the measurement function. Many aspects in this 
discussion are completely analogous to the 3-jet case, so 
we will focus on those where the extension to more jets is 
nontrivial, which are the kinematic dependence and the 
color structure. 



A. Kinematics and Power Counting 

For our observable we again use 7V-jettiness defined in 
Eq. (|4.ip . The two beams are included using two refer- 
ence momenta qa and qh, which correspond to the mo- 
menta of the incoming partons, and the corresponding 
dimensionless reference vectors qa,b — <la.b/Qa,b, which 
determined the separation between the beam and jet re- 
gions (see Refs. 0,0] for more details). We will consider 
the cross section differential in the iV-jettiness contri- 
butions 7^, 7fc, 7i, . . . , T/v, where Ta.t measure the con- 
tribution from the beam regions. In addition we mea- 
sure the small dijet invariant mass t and the energy ratio 
z — El /{El + E2) for the two nearby jets. 

As before, we label jets 1 and 2 as the two nearby jets, 
and consider the limit in Eq. (|2.12|) , with all jet energies 
parametrically of the same size, such that we have 

t = S12 < Sij ^ , st = §12 < s,y 1 . (6.1) 

(To have a manageable notation, we specify 1 and 2 to 
be two final-state jets. The case where jet 1 is close to a 
beam, such that sia = 251 • qa ^ Sy is completely anal- 
ogous and does not involve additional complications.) 

The power expansion in = t/Q^ is again defined by 
choosing a common reference vector qt for jets 1 and 2, 
as in Eq. (14. 8p . This gives 

SQ, = 2qt ■ = + OiXt)] = S2,[l + 0{Xt)] (6.2) 

for each jet i 1,2, which generalizes sq from the 3- 
jet case. The corresponding dijet invariant masses in the 
SCET above ^/t are then given by 

Sli — QlQi SQi , S2i = Q2Qi SQi , 

SU = {SU + S2^) = [{qi + q2) + q^?[l + 0{\t)] , (6.3) 



which we can also express as 

su^xsti, S2i = [l - x) sti , = — -^i— . (6.4) 

Wi + W2 

For the geometric measure, Qi — Ei, we have x = z as 
before. 

The factorization of the measurement function fol- 
lows the same logic as in the 3-jet case. By tak- 
ing the collinear, csoft, and usoft limits of the full A'^- 
jettiness measurement function, we obtain the general- 
izations of the measurement functions in the 3-jet case. 
The collinear and csoft measurement functions are not 
affected by the presence of additional jets and so are un- 
changed from the 3-jet case. The generalization of the 
usoft measurement function for the iV-jet case is given 
by 

Mi^iTr,Tr,Tr,---,T;}n (6.5) 

= n ^i^^' ~ • n ^c^"' - . pr) , 

i=l,2 i#l:2 

where the momentum operator P"'* is defined in 
Eq. (|4.30p . and the 6"''(p) are defined by the obvious 
generalization of Eq. (I4.23|) . 

&rip)= n m-p-qfP)o[cosMp)], 

0r(p)- II 9{q,-p-qfp)e[-C0sMp)]: (6.6) 
j#l,2 

Qr^lAP)=^i^t■p~q^■p) n e{q,-p-q,-p). 

The reference vectors qi for i 7^ 1 , 2 can now have a 
nonzero component in the transverse plane, and can 
therefore have 4>t dependence. This implies that the jet 
regions for jets 1 and 2 are in general not symmetric 
(unlike the 3-jet case, where they were symmetric up to 
power corrections in t/Q"^). 

B. Factorization 

The hard factorization in SCET_|_ proceeds through the 
same basic steps as for the 3-jet case. We first match from 
QCD to SCET at the hard scale Q, 

Xqcd(2 ^ ^-1) = (7V-l|djv-i(M)|2)Cjv-i({s.,},Ai) , 

(6.7) 

where A^qcd(2 ^ iV - 1) is the 2 ^ - 1 QCD am- 
plitude for the process we are interested in, and On-i 
is the corresponding 2 ^ N — 1-jet operator in SCET, 
discussed below. We will again use k to denote the de- 
pendence on a specific partonic channel when needed, 
but to simplify the notation, we mostly suppress the la- 
bel K in what follows. As before, the bare loop diagrams 
of Om-i vanish in pure dimensional regularization, so in- 
cluding counterterms the renormalized matrix element of 
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On-i equals the tree-level result plus pure 1/e IR diver- 
gences which precisely cancel against the IR divergences 
in the QCD amplitude. Therefore, to all orders in per- 
turbation theory, Cm-i is given by the finite parts of 
7Wqcd(2^^-1). 

The operator O^-i in the matching in Eq. (|6.7p has 
the form 

C'jv-l ^ ^ ['^"t] [^"a] ['^"b] [^"3] ' ' ' ['^njv] 

X [y„,r„^r„,K„3---r„„] . (6.8) 

We let C„. denote a (usoft-decoupled) gauge-invariant 
coUinear field in the Ui direction, which can be a coUinear 
quark, antiquark, or gluon, and F represents the spin 
structure connecting the different fields together. In gen- 
eral there are many such structures possible, so Eq. (16. 8p 
really represents a set of operators. As before, jets 1 and 2 
are described by a single coUinear field C„j in the rit direc- 
tion, and are the fields for the incoming partons, 
and C„3 to C„„ are the fields for the outgoing partons 
that initiate the remaining final-state jets for i > 3. The 
usoft Wilson lines are written generically as Yn^ without 
any reference to their particular color representation. 

The operator On-i and Wilson coefficient Cn-i in 
Eqs. (|6.7I) and (|6.8p are now vectors in the color space 
spanned by the 2 -I- TV — 1 external partons, as indicated 
by the vector symbols. That is. 



N- 



— r)t"i 

-1 = 



1 



(6.9) 



where ai is the color index of the ith external particle. 
The product of all usoft Wilson lines in Eq. (16. 8p is a 
matrix in the same color space, 



1 0t ■■■/3jv lof-tiN 



(6.10) 



The vertical bar separates the column indices (on the 
left) and row indices (on the right) of the matrix. The 
color charges act in the external color space as 



(T^C)--^'-- 



(6.11) 



where the three lines are for the kth particle being an out- 
going quark or incoming antiquark, an incoming quark or 
outgoing antiquark, or a gluon, respectively. The prod- 
ucts Tj • Tj ~ T^T^ are matrices in color space. 
From Eq. (|6.1ip it is clear that Tj for different i com- 
mute. 

In the next step we match from SCET to SCET+ at the 
scale \/t. From the construction of the effective theory 
in Sec. IIIIl it should be clear that the relevant operator 
in SCET+ is constructed out of -I- 2 coUinear fields, for 
the two incoming and N outgoing partons in the hard 
interaction, csoft fields that interact with the coUinear 
fields in directions 1 and 2, and usoft fields that interact 



with all coUinear degrees of freedom. The 2 
operator in SCET+, O^, is obtained from Eq. 
the analogous replacement as in Eq. p.33p . 



[Cni]^ [Cna]^ [A„j^ A"„2 Tt 



|/3i/32|at 



3.12) 



In the second line we performed the csoft decoupling 
Eq. (|3.25|) . which produces the csoft Wilson lines Xm 
and (dropping the superscripts that distinguish the 
fields before and after the field redefinition). The color 
generator Tj is contracted with the color indices of the 
daughter fields as shown in the first line of Eq. (|6.12l) . 
From the product of csoft Wilson lines we define 



* -|_^a---/3iv |a„---Qjv 



(6.13) 



which is a color space matrix that takes us from 
{N + l)-parton to {N + 2)-parton color space, and 
1,3„---/3„|q„---q„ ^ ^/3„a„ . . . j/3«a« ^enotcs the iden- 
tity in the remaining iV-parton color space for partons 



a, 5, 3 
6 



N 



. , N . The operator OJj then has the form 

= F [C„J [C„,] [C„J • • • [C„„] [X] [Y] , (6.14) 



where the product of all coUinear fields is a row vector 
in {N + 2)-parton color space. We have grouped the 
different factors in square brackets belonging to different 
sectors, which do not interact with one another through 
the leading-order SCET+ Lagrangian in Sec. IIIIl 

The matching from SCET to SCET+ takes the form 

{N\6i_,(pp) = {N\d+Hfi)\2)C+{t,x,fi). (6.15) 

Since On-i and are both vectors in color space, in 
principle C+ could be a matrix in color space. However, 
since the different sectors in both SCET and SCET_|_ are 
explicitly decoupled, the matching coefficient C+ is actu- 
ally determined by the 1 — > 2 matching in Eq. (|6.12p 



^7it (i , .x) C7ij^C7i2 A^ij^ V^i^ 



(6.16) 



In other words, C+ = is universal and only depends 
on the specific 1 — >■ 2 splitting channel q ^ qg, g ^ gg, 
or g qq. In Appendix |Al we use reparameterization 
invariance to show that C+ depends only on t, x, and the 
azimuthal angle of the splitting. 

Using the same arguments as in Sec. IV A] we can relate 
C+ to the coUinear limit of the 2 — >■ A QCD amplitude. 
Since the matching onto Ojv-i in Eq. (|6.7p is independent 
of the external state, we have 



A^qcd(2-> A^) 
= {N\dl_,\2)C, 



{N\d 



N-1 



+^|2)Cjv-i({^ 



,})C+(i,x), 



(6.17) 



where A^qcd(2 A^)|t<ss,j is the 2 ^ A^ QCD ampli- 
tude expanded in the coUinear limit of partons 1 and 2 be- 
coming close, and in the second step we used Eq. (j6.15l) . 
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The loop diagrams of again vanish in pure dimen- 
sional regularization, so the product Cn-iC+ is given by 
the IR-finite part of A^qcd(2 — )■ N)\t<g:sij ■ On the other 
hand, as we saw above, Cjv-i contains the IR-finite parts 
of A^qcd(2 -> iV - 1). Therefore, 



M'^'^ui'^ ^N) ^ C+{t, x) Xjj^cD "(2 ^ - 1) . 

(6.18) 

It is well-known that the iV-point QCD amplitudes in the 
coUinear limit t <^ Sij factorize into (N — l)-]X)int ampli- 
tudes times universal splitting amplitudes j36l44ll |. Just 
like the IR-finite parts of the full amplitude determine 
the hard matching coefhcient, it follows from Eq. (|6.18p 
that the same is true for the splitting amplitudes: The 
IR-finite parts of the splitting amplitudes directly deter- 
mine the matching coefficients C+ for the different par- 
tonic channels. Taking the square of the one-loop results 
for the q ^ qg splitting amplitudes from Ref. [38] and 
summing over helicities reproduces the expression for 



in Eq. (|5.1ip . In the same way, the one- loop results for 
C+ for the other splitting channels can be obtained. 

The cross section in SCET-|_ is obtained from the for- 
ward matrix element of in Eq. (|6.14p with the mea- 
surement function inserted, 



da^\C+\'Cl_,{2\Y^X^l[cl 

i 

X MN{{Tk}) '[lCn^XY\2)CN- 



(6.19) 



Using the factorization of the measurement together with 
that of the operator, the matrix element factorizes into 
independent collinear, csoft, and usoft matrix elements. 
The collinear matrix elements produce N jet functions 
and two beam functions, which are all diagonal in color 
and contribute a factor of 1 = Yii S"^^'- . The remaining 
soft matrix element is given by 



Qi^Pt-PN ^Q^Y^ l3t---0N\l3't---l3'N xf l3't---l3'N\a"a'2---a'N J\^'U^ (6.20) 

where we explicitly wrote out the color indices in the product of csoft and usoft Wilson lines. From Eq. (j6.13l) we 
know that X is diagonal in color except for the 1,2,< subspace, so the product X^'X has only two nontrivial color 
indices P't\a[. The only object we can form from these is J^*"*, which implies that the csoft matrix element is entirely 
color diagonal, 

{0\[x'<M''%ki,k2) Xf''"'^"^"'"'""' \0) = S'+(fci,fc2)l'^*-^"l"*-"" . (6.21) 
The csoft function S-^- is the same as in the 3-jet case, 

s^{h,k2,fi) = ^tr(o|r[K!,Tjx4xtjx-(fci,fc2)r[x„,x„,TtKj|o), (6.22) 

where we restored the proper time-ordering, the trace is over color indices, and the color normalization constant, c", 
is such that at tree level S'^(ki,k2) = 5{ki) 5{k2)- Like C^, the csoft function is universal and only depends on 
the color representations of the partons l,2,t involved in the 1—^2 splitting. The explicit form for qt — > qig2 was 
given in Eq. (|4.36p . Using Eq. (|6.2ip in Eq. (|6.20p . the remaining usoft matrix element yields the usoft function 

S%_^{{h},^i) = -^{Q\f[Y^]Ml^[{h})T[Y]\0) , (6.23) 
which is a matrix in + 1-parton color space, and the color normalization factor, c^_x, is such that at tree level 

%-i({fcO) = ia'5(fco- 

Having discussed the color structure, assembling the full factorization theorem for the iV-jet case now follows the 
usual steps. For the cross section differential in the 7i, the small dijet invariant mass t, and the energy fraction z, we 
find 

da 

ATa d7b dTi • • • ATn dt dz 

El 



E1+E2 



jd\d^L{q) jd^N{{q^}) Mn{^n. $l) {2Tif5^(qa + qb-Y,q.- q) S{t - s^) S(^z 

i 

X Jdxadxb JdsadsbBf^^{sa,Xa,n)B,,^{sb,Xb,n)Y]_ jds^Jf,Xsi,tJ^)\C1{t,z,^i)\'^jdki dk2S1{ki,k2,n) 

K i 

X C;^li($7V, $l,m) S^_i (Ti - ^ - kuT2 - ^ - k2,Ta - . ■ ,Tn - M)^^^-i($Ar, $i,/x) . (6.24) 
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Here we have included the possibihty for the hadronic 
system to recoil against a color-singlet final state with 
total momentum q and internal phase space 3>l (9) ■ This 
allows us, for example, to describe W/Z + jets at the 
LHC. The massless iV-jet phase space for the N final- 
state jets is denoted as $Ar({gi}). The incoming mo- 
menta qafi are given by qa^b = a;a,&£'cm(l, ±^)/2, where 
Ec-m is the total hadronic center-of-mass energy, z points 
along the beam axis, and are the light-cone momen- 
tum fractions of the colliding hard partons. The restric- 
tion Mm enforces any phase space cuts on the final state, 
such as requiring that the jets be energetic and only one 
dijet invariant mass be small. The sum over n again runs 
over all relevant partonic channels. The jet functions, 
J^i , and beam functions, i?^^ ^ , arise from the collinear 
matrix elements as described above. The beam functions 
depend on the momentum fractions and describe the 
collinear initial-state radiation of the incoming hard par- 
tons |3t42|, whereas the jet functions describe collinear 
final-state radiation of the outgoing hard partons and are 
the same as in the 3-jet case. 

Since each of the ingredients in Eq. ()6.24p only depends 
on a single scale, the standard RG procedure can now 
be followed to resum the large logarithms in the cross 
section, by evaluating each function at its natural scale 
and then evolving all functions to a common scale /i. 
The solutions to the RG equations are straightforward 
generalizations of the explicit results presented for the 
case of e+e~ 3 jets in Sec. IVIII 



C. Consistency 



and given by the IR- finite terms of A^qcd(2 — >■ N)^ since 
as before the loop corrections to the bare matrix element 
of Ojv vanish in pure dimensional regularization. Thus, 
from Eq. (|6.17l) and the discussion below it, it follows 
that 



CN{{sij},^i) = Tt CN-l{{siJ},^l) C+{t,x,fj,) . 

(6.27) 

Note that Cn and Cn^i are vectors in the color space of 
2 -I- iV and 2 -|- — 1 external partons, respectively. We 
can write this with explicit color indices as 

(^aia....a„ ^ rj.^l^2\a, ^a.^-a^ C+{t,X,ll) . (6.28) 



Taking the derivative with respect /i and using 



(6.29) 



we get 



fi—C+{t, X, ^) = 7C+ (i, X, ^) C+{t, X, pi) , 



■Tl7c„(K},/i) T, 
7Civ-i({sii},M), 



(6.30) 



Having derived the factorization theorem for pp N 
jets, we now discuss its renormalization group consis- 
tency. We already know the structure of the anomalous 
dimensions of the new objects C+ and S+ from Sec. IV Dl 
Nevertheless, it is instructive to see how they fit together 
with the anomalous dimensions of Cn-i and Sn-i in- 
cluding the more complicated color structure now. 

We will follow the same logic in our discussion here as 
in Sec. IV Dl Our starting point is the factorized pp — > A^- 
jet cross section for A^-jettincss in SCET given in Rcf. j7|], 
which has the form 



E 



N 



® [C^^S^Cl,] . (6.25) 



We then argue that the combined anomalous dimensions 
of C+Cn-i and S+Sn-i must reproduce those of Cn 
and Sn in the limit t <^ Sij. Using this, we rederive the 
general form of the anomalous dimensions for C+ and 
5+, reproducing those found in Sec. IV Dl 
1. Consistency of the Hard Functions to Determine 'yc+ 



The Wilson coefficient Cn hi Eq. ()6.25p is obtained 
from 



where T^Tj — ICt, with Ct determined by the parent 
parton splitting to 1 and 2, Cq = Cq = Cp and Cg = Ca- 
Since C+ is a scalar function and does not know about 
the color space, the difference between the anomalous 
dimensions in the second equation must be proportional 
to the color identity. 

The /^-dependence of ■^Cn (m) has the all-order struc- 
ture m, M, S-m 



7Civ(Ai) = -rcusp[a3(M)] E \ ^ 



X In 



^~^^^"r^"''° +7c.[a.(M)], (6.31) 



where A^- = 1 if both i and j are incoming or outgo- 
ing and Aij — otherwise; Sij is always positive. We 
do not make additional assumptions about the all-order 
structure of the noncusp term 7c„ [as] i which can in gen- 
eral be a matrix in color space. Beyond two loops and 
for four or more external directions it can also depend 
on the Sij through conformal cross ratios of the form 



Mqcd{2^N) (iV|0]v(/i)|2)Cjv({s.,},/i), (6.26) SijSki/ SikSji [4J,|45[. Hence, we can write 7Cn-i as 
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7Civ-i(M) = -rcusp[as(Ai)]<^ 2^ Tf • T,ln ^ + — 2 [ + 7Civ~i [a^W] ■ 

(6.32) 

Using Eqs. (|6.3p and (|6.4p to expand in the limit i ^ Sy , we find 

+ i (t? - - Tl) In - i (t? + T? - T^) In X - i (t? - + T^) ln(l - x) | 

+ 7tTl7c„K(M)]Tt, (6.33) 

where we have used the identity J2i^i 2 = ^(Ti + T2) and the following relations 

-^Tj(Ti+T2)-T,Tt =TfT,, ^^tJ Ti • T2 = i (T^ - T? - T^) , (6.34) 

to write the contribution of 7(7„ entirely in 2 + — 1 parton color space. Equation (I6.34p generalizes Eq. (I5.44p to 
the A^-jet case. Taking the difference of Eqs. (|6.33p and (|6.32p . we obtain the result for jc+ valid to all orders in 
perturbation theory, 

7C+ (t, x, /x) 1 = -rcusp[a.(M)]^ (t? ~Tj- T^) In ^^-^ + 7^^ [a,(^)] 1 , 

7c+[as] 1 = -r,usp[as]|-^(T? + - T^) Inx - ^(t^ - + T^) ln(l - : 

+ (7c[«.] + 7c[«.] - I'M) 1 . (6.35) 

Here we have used that the difference of the noncusp terms in the coUinear limit must be diagonal in color to all 
orders since as above in Eq. (j6.30p jc+ [os] must be color diagonal, 



^ Tl7c„[a,]Tt -7c„„Ja.] = (7^^] + 7^^] - 7^[«,]) 1 . (6.36) 



This condition provides a nontrivial constraint on the collinear limit of 'JCn ^^id is equivalent to what was used in 
Refs. [13, SI] in deriving constraints on the form of jcn beyond two loops. (It cannot be spoiled by conformal cross 
ratios appearing in 'JCn > which would require at least four distinct momenta and thus cannot be reduced to functions 
of t and X only.) 

To see that Eq. (j6.35p agrees with the result for ^0+ in Eq. (j5.45p . we note that in the three-parton color space used 
in Eq. (I5.45P we have T3 = — (Ti + T2). Equivalently, we can write the color factors in Eq. (|6.35p in the external 
+ 2 color space using the inverse of Eq. (I6.34p , which amounts to the replacements 

^^T2-T2~t2) ^Ti-T2, ifx^ + T^-T^) ^Ti-(Ti+T2), ^ (t^ - + T^) ^ T2 • (Ti + T2) . (6.37) 



With these replacements, Eq. (j6.35p agrees with the anomalous dimensions of the splitting amplitudes given in 
Ref. [33] . Note that in iV + 2 parton color space is not color diagonal but essentially proportional to T(TJ. 
Working in the N + 1 parton color space has the advantage that 7c+ is manifestly color diagonal as in Eq. ()6.35p . 



2. Consistency of the Soft Functions to Determine 75^ 

The cross section for pp ^ N jets in SCET_|_ must reproduce the corresponding cross section in SCET expanded in 
the limit t <^ sij. Since the product of hard matching coefficients in SCET+, C+Cn-i reproduces Cn in SCET, as 
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we saw in the previous section, and the beam and jet functions are the same in both theories, the soft functions must 
satisfy 



SNi{ki},^i) 



^ J dk[dk2 S+{k[,k2,^i) Tt S^-iiki - k[,k2 - fcj, fc^, . . . , fcjv, A^) tJ , 



(6.38) 



which is the iV-jet generahzation of Eq. (|5.47p . Taking the derivative with respect to fi and using the known evolution 
equations for Sn in SCET from Ref. 0, we find 

li-^S+{ki,k2,n) = J dk[dk'2^s+iki - k[,k2 - k'2, fi) S+{k[,k'2, fJ.) , (6.39) 



is. 



+ (/ci,fc2,Ai)i n ^ik^)^^^^ll[lsA{k^},^l)+%J{h},^i)] 

«#1,2 * 



The Tj7SjvTt projects the anomalous dimension of Sn on the 2 + — 1 parton color space. Similar to what we saw 
for C+, the difference on the right-hand side of the second equation must be color diagonal, since knows nothing 
about the larger color space. 

We can now use the all orders form of the A^-jettiness soft anomalous dimension derived in Ref. Q, 



lsA{k^}, ^i) = -2rcusp[as(M)] • 



Co(^^)~'-^A,,S{h) Y[6{k„.) + jsMf^)]l[5{k,). 



to determine in the limit St <C s^j, which corresponds to setting su — S2i — sq^ [see Eq. 
straightforward to extract the anomalous dimension of 5*+ using Eq. (|6.39p . We obtain 



(6.40) 
It is then 



ls+ikuk2,^^)l = -2rcusp[«.(M)]^(T? -T2^ ^ 



Co 



J^)sik2)+Co(^)sik,) 



k2 



1S+ K, kiM 1 = -re„,p[a,](T2 - T?) 



Co 



(|)5(fc.)-/:o(|)<5(fci) 



(7}[a.] + 27i,[a,] + 7,2 [a,] + 27^ - 7^[«,] - 27*,[«,]) 1 5{k^) ^(fca) 



+ 7S+ [Q!s(/z),fci,fc2] 1 , 



3.41) 



This result agrees with Eq. (|5.58p for e+e~ — >■ 3 jets as it must since the csoft function is universal for a given splitting 
channel. 

To arrive at Eq. (|6.4ip . we have used that the difference of the noncusp terms in 7s„ and 7Siv-i in the collinear 
limit must be color diagonal because 75,^ is, i.e.. 



^t\^sM Tt - 7S„_Ja.] = -{l)Ws\ + HcWs] + 7,7[«s] + 27^ - 7* - 27*,[a,]) 1 . 



3.42) 



This condition is equivalent to Eq. (j6.36l) . since 75'„ and ICn related by consistency. We have also used the 
relations in Eq. (|6.34p to rewrite the different color structures appearing in 75'„ . The terms proportional to Ti • T2 
in 75„ directly give the terms proportional to (Tj — — T2)/2 in the first line of Eq. (|6.4ip . The terms in 75^ 
involving neither Ti nor T2 immediately cancel with the corresponding ones in Isn-i- The nontrivial terms in 75^ 
are those proportional to Ti • and T2 • with i 7^ 1, 2, which we can rewrite as 



Ti -T, 



1 



T2 • T,: 



1 



Co(^^)s{h)+Co(^^)s{ki) T\ Sikj) 



(Ti+T2)-T, 



k2 

1.. 



S{k,)+£o 
ki 



h 

ki 



SQ^^J■ 

Ti + T2) • T, -^^Co(^^) n '5(%) 



5ik2) + 



Sik2) 

k2 



n -^Ifc.) 



Siki 



n '^(^^■) 
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i(Ti-T2)-T,i 



(6.43) 



r 



After projecting (Ti+T2)-Ti TfTi, the first two lines 
on the right-hand side cancel against jsn-i- To see this 
first note that as in the 3-jet case, our SN-i{ki, k2, ■ ■ ■) 
is related to the usual SN-i{kt, ■ ■ ■) from Ref. by pro- 
jecting onto kt = ki + k2. To determine how the fcj 
measurement is split between jets 1 and 2 notice that 
for general N, SN-i{ki,k2, ■ ■ ■) is no longer symmetric 
under fci o fe, because the boundary to the ith jet re- 
gion can be different for jet regions 1 and 2. However, 
the fact that the anomalous dimension of S-^- must be 
color diagonal requires that the first line on the right- 
hand side cancels against JSn-i- This tells us that the 
contribution to the anomalous dimension of SN-i{kt, ■ ■ ■) 
involving Tj • TiCo{kt) must again be split up symmet- 
rically into Tt ■ T,[Co{ki)Sik2)+Coik2)Siki)]/2 in 7s„_, 
in our case, just as for 75^ in Eq. (j5.56p .'^ Finally, since 
the last line on the right-hand side of Eq. (|6.43l) has no 
dependence on i other than itself, we can sum over 



-(Ti + T2). This yields the term 



proportional to - in Eq. ([OT|) . 



VII. NUMERICAL RESULTS 

In this section we present some numerical results for 
our example of e^e~ — > 3 jets. For simplicity, we project 
onto the total 3-jettiness of the event [see Eq. (14.31) ]. 

r = rN=3^ri + r2 + r3. (7.1) 

Just as thrust characterizes how 2-jet-like an event is, the 
total A/'-jettiness of an event characterizes how A^-jet-like 
an event is, and can be used as a veto against additional 
jets. The factorized cross section for T is obtained from 
Eq. (|4.35p by projecting onto T using Eq. (|7.ip . 



da 



dT dt dz 



(7.2) 



= S E ^2(g', ^i)Hl{t, z,^i)l[^ds, J,, {s,,^i) 

K i 

x/d.54<..M)&(r-|L_g_|L_.,, 



To see that 7s„_j must still be symmetric under fci <-> k2, we 
can split up the phase space region associated to jets 1 and 2 into 
a </<i (p) independent region enclosing qi and (j2 and the UV-finite 
difference to the boundaries of the ith jet region. In Ref. [7| it 
was demonstrated at one loop that the UV divergences associated 
with jet regions i and j only arise from i,j hemispheres, whereas 
the contributions depending on the boundary with the remain- 
ing jet directions are UV finite and do not affect the anomalous 
dimension. 



The different hard, jet, and soft functions were discussed 
in detail in Sec.|Vl They are renormalized objects, and in 
the cross section must all be evaluated at the same scale 
fi. To resum the large logarithms in Eq. (j7.2p we have 
to evaluate each function at a scale where its perturba- 
tive series does not involve large logarithms and then RG 
evolve all functions to the common scale /z, 



da 



dT dt dz 

- ^ 
^ n2 



X H'lit, z, ^lH+) W^^it, z, ^i) 

^11/ 'i^^ds ■ J«. (sj - s ■, fij) Uj^^ (s ■, fij, ^) 

X / dkdk' S+{k- k',fis+)Us+{k',^is+,fJ-) 



X ld£S2(T-^-^-^~k-£,fis2 

Wl W2 W3 



X f7s2(^,MS2;M) 



(7.3) 



The evolution kernels Ux for each function in Eq. (|7.3p 
are given in Appendix [C] After RG evolving the hard, 
jet, and soft functions, we obtain a distribution with all 
large logarithms of T and rrijj resummed. 

We are interested in the distribution in rrijj , for which 
we integrate the cross section over T and z, 



da 

dm 11 



2m 1 



dz 



dT 



da 



dT dt dz 



(7.4) 



where the factor of 2mjj is the Jacobian from changing 
variables from t = m?,- to rrijj, and we suppress the de- 



pendence of da/dmjj on Tut and Zcut- For our numerical 
analysis, we choose 



Q = 500 GeV , 7;ut = 10 GeV , 



(7.5) 



The Qi are given in terms of z and t in Eq. (|4.17p . For 
this Zcut, the ratio of nearby jet energies ranges between 
0.5 and 2. For T^ut = 10 GeV the two nearby jets have 
a typical jet mass of 35 GeV. In terms of rrijj and T^ut, 
the scales for the different functions in Eq. (|7.3p have the 
scaling 



Tut 
y/st 



MS2 



Tcu 



\/ QiTx_ 



(7.6) 



The simultaneous resummation of logarithms of Tut 
and rrijj is challenging because the hierarchy between 
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the scales depends on several kinematic variables. At 
large rrijj we are in the kinematic region of 3 equally 
separated jets. In this 3-jet limit, power corrections 
in m'jj/Q^ become important and the resummation in 
rrijj must be turned off, which requires iih+ — M_ff2 ^^'^ 
fis^ — ^J■S2■ For any mjj we also have the requirement 
that fJ.H+fJ'S+ ~ MJiMJ2- To satisfy these requirements 
we choose the scales as follows: 









= \/ QiTcnt 















'IJ-S2 



fJ^runiQ, rrijj) 



(7.7) 



where sq = Q^/lQaiQi + Q2)]- The scales ^h, lJ.j,, 
and 11S2 do not depend on nijj so we simply use their 
canonical scales in Eq. (|7.6p . The -dependent scales 
^^H+ (iTijj ) and 113+ {mjj ) are given in terms of the profile 
function 



2a rrijj + fiQ 

(rrijj — 7712)^ 
fi — a- 



m2 — mi 



rrijj < mi , 

mi < mjj < 7712 , 



mjj > 7712 



Mo 



7771 + m2 



(7.8) 



such that they have a quadratic approach to and 
over the region mi < mjj < 7772. Analogous profile func- 
tions have been used previously in different contexts [46l - 
lisj . We choose the parameters of /ij-un as 

7771 = 200 GeV , 7772 = 300 GeV , = GeV . (7.9) 

To get a measure of how important the resummation in 
ma is, we will compare to results with the resummation 
turned off. This is achieved by setting 
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in 777, 



+ 



2 1 



Ms 



+ 



MS: 



2 ' 



(7.10) 

with Hh2^ MJiJ MS2 as in Eq. (|7.7p . Note that the distri- 
bution without resummation in mjj is what one would 
obtain from directly using SCET with 3 jet directions 
(up to the power corrections in t/Q'^ which we have not 
included). 

The approach to the 2-jet region where the two nearby 
jets merge into a single large jet is more complicated. 
This happens for mjj < mj ^ 35 GeV, i.e. when the 
dijet invariant mass becomes smaller than the individ- 
ual jet masses. Here, ^.r^ and become equal to 
and eventually becomes larger than ■ The 3-jet 
observables 777^^ and 3-jettiness are not meaningful any- 
more once the two close jets merge into each other. In 
addition, the proper theory would now be SCET with 
two collinear sectors. We leave a more detailed investi- 
gation of this transition to future work. For illustrative 
purposes, we will plot our results all the way down to 
mjj = 5 GeV, where becomes equal to ^g^. 
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FIG. 5: The distribution in mjj for Q = 500 GeV, T^ut = 
10 GeV, and Zcut = 1/3 as in Eq. 1)7. 5|) . The bands show the 
perturbative scale uncertainties at NLL (light blue) and NLL' 
(dark orange) as explained in the text, with the central values 
given by the center of the bands. The dots show the result 
obtained from Pythia. The dotted line indicates the value 
rrijj — 35 GeV below which we enter the 2-jet region where 
our expansion breaks down. 



A detailed numerical study is beyond the scope of this 
work. Below, we will present results at next-to-leading 
logarithmic (NLL) order, which uses the one-loop non- 
cusp and two-loop cusp anomalous dimensions in the run- 
ning and tree-level matching, and at NLL' order, which 
combines the NLL running with the one-loop matching 
corrections. It is straightforward to extend our results 
to next-to-next-to-leading logarithmic (NNLL) order. To 
obtain a rough estimate of the perturbative uncertainties 
we vary /iff2 7 MJiJ MS2 by factors of two. By varying ^h2 
and ^S2 automatically let iih+ and vary accord- 
ingly. This tends to give the largest effect in the variation 
of /i-ff+ and j whereas varying the parameters in ^^un 
within reasonable ranges has smaller effects. We then 
take the envelope from varying each of fj,H2 j M Ji i M J2 1 
/ij.,, /is2 individually while keeping the other scales fixed 
and in addition from varying all of them up and down at 
the same time. The largest variation mostly comes from 
the individual variations of /igj and the gluon jet scale. 

In Fig. [5l we plot the NLL and NLL' distributions from 
Eq. (|7.4p for our default scale choices. The distribution 
is stable in going from NLL to NLL', with the expected 
reduction in the scale uncertainties. The ninja region, 

cor- 



which has the hierarchy of scales shown in Fig. |l(b) 



responds to the region to the right of the dotted line in 
the plot. For smaller mjj we enter the 2-jet region where 
our expansion breaks down. For mjj ^ 150 GeV we enter 
the 3-jet region where power corrections of order mjj/Q 
become important. 

In Fig. [5] we also compare our resummed predictions to 
the distribution obtained from PythiaS [i^ . To make the 
distribution in Pythia, we simulated e'^e~ — )■ qq events 
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FIG. 6: Left panel: The distribution in rrijj for Q = 500 GeV, Tcut = 10 GeV, and Zcut = 1/3 as in Eq. (|7.5|) . Right panel: 
Same as the left plot but shown as the percent difference relative to the NLL' central value. The bands show the perturbative 
scale uncertainties as explained in the text. The dark orange band shows our result including the resummation (resum.) in 
rrijj at NLL'. The medium blue and light green bands show the results at NLL and NLL' with the resummation only in Tcut 
but not in m,, . 



and determined the reference momenta qi from the first 
hard emission. For simplicity we use this as an approx- 
imate alternative instead of determining the qi by fully 
minimizing the 3-jettiness for each event. (We found that 
using a jet algorithm to determine the reference momenta 
introduces a significant bias as the nearby jets get close 
to each other and the events become more 2-jet-like. In 
this limit, jet algorithms tend to merge nearby clusters 
of high energy into a single jet [50| . and hence become 
unsuitable to determine the qi.) Pythia resums the kine- 
matic logarithms of rrijj in the parton shower, which is 
formally correct at leading-logarithmic order, but also in- 
cludes various other effects that contribute at NLL. We 
see that the Pythia distribution agrees well with both the 
NLL and NLL' distributions even down to rrijj ~ 25 GeV. 
Note that Pythia is producing exclusive samples of events 
and does not report any systematic uncertainties. This 
makes it challenging to interpret the theory uncertain- 
ties in the distributions generated by the Monte Carlo. 
Even though the NLL and NLL' distributions are more 
accurate than the parton shower, since Pythia contains 
many NLL effects it is reasonable to take the uncertainty 
band of the NLL curve as a proxy for the systematic un- 
certainties in the Pythia distribution in the range of rrijj 
where we trust our calculation. 



In the left panel of Fig. [SJ we compare the NLL' dis- 
tribution including the resummation in rrijj to the NLL' 
and NLL distributions without rrijj resummation. In the 
right panel we show the same results in terms of the per- 
cent difference from the central NLL' curve with rrijj 
resummation. We see that for rrijj < 75 GeV the re- 
summation in rrijj becomes important, where the results 
without rrijj resummation become unstable due to the 
presence of unresummed \n{mjj /Q). 



VIII. CONCLUSIONS 

In this paper we have constructed a new effective field 
theory, SCET4-, which can describe exclusive multijet 
events at the LHC where the kinematic configuration of 
final-state jets give rise to a hierarchy of scales. We fo- 
cused on the case where the dijet invariant mass rrijj of a 
pair of jets is much smaller than all other dijet invariant 
masses, which are of order the total center-of-mass energy 
Q. This results in the presence of three separated scales 
m <C rrijj <C Q, where m is the typical jet mass. We have 
shown that using SCET+ we can simultaneously and sys- 
tematically resum the logarithms of both the individual 
jet mass, ln(m/(5), as well as the dijet mass, ln(mjj/Q), 
to higher orders in perturbation theory, which has a wide 
range of applications. 

Separating the additional kinematic scales that arise 
requires a successive hard matching from QCD onto 
SCET with TV - 1 jet directions and then onto SCET+ 
with N jet directions. The matching from SCET onto 
SCET+ introduces a new hard matching coefficient, 
which encodes the splitting to two nearby jets, is univer- 
sal for a given splitting channel, and depends only on the 
scale rrijj to all orders in perturbation theory. We showed 
that it is given by the finite parts of the known universal 
splitting amplitudes of QCD in the collinear limit. 

SCET+ is an extension of SCET with an additional 
csoft mode that has virtuality m'^ /rrijj and is necessary to 
describe the soft radiation between the nearby jets. The 
usoft modes of SCET have virtuality m'^/Q^ and there- 
fore cannot resolve the two nearby jets. Thus, SCET+ is 
required to properly separate the physics between these 
two scales. We constructed the Lagrangian of SCET+ 
and showed how to decouple the collinear, csoft, and 
usoft interactions. We discussed the structure of gauge 
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invariance in SCET_|_ and used this to construct gauge in- 
variant iV-jet operators, which require new csoft Wilson 
hues. The factorization of the csoft and usoft interactions 
leads to the factorization of the regular soft function into 
a new universal csoft function, S+, and a usoft function. 

As an example of our new effective theory, we derived 
the factorization for the dijet invariant mass spectrum 
of e+e~ — >■ 3 jets, using A^-jettiness to define the jets. 
We determined all contributions to the factorized cross 
section at NLO and derived the form of the anomalous 
dimensions of C+ and 5*+ to all orders in perturbation 
theory. We also extended our factorization theorem to 
pp ^ N jets plus leptons, discussing in detail the kine- 
matic scales that appear and the nontrivial color struc- 
ture for N jets. We demonstrated that the new ingredi- 
ents, C+ and S+, are color diagonal and independent of 
the number of additional jets in the process. 

We gave explicit numerical results for the example of 
e+e" — 3 jets using 3-jettiness, T, resumming both the 
jet-mass logarithms of T/ Q and the kinematic logarithms 
of rrijj/Q to NLL'. We found that the resummation of 
rrijj is important over a wide range. Our results show 
good agreement with Pythia, which resums the kinematic 
logarithms of rrijj at leading order. 

SCET+ has a wide range of applications, and expands 
our ability to make accurate theoretical predictions rele- 
vant for jet-based new-physics searches at the LHC. This 
includes cases where dijet invariant masses lead to a nat- 
ural scale hierarchy, ranging from the study of jet sub- 
structure to improving parton-shower programs, which 
are currently the only tool available to resum logarithms 
of kinematic scales in exclusive jet cross sections. 
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Appendix A: All-Order Scale Sensitivity of the Hard 
Function 



We have argued that the matching coefficient 
C+{t, z, /i), which arises from matching SCET to SCET+, 
can only depend on the scale t, which is the scale where 
the matching is performed. In Sec. IV Al we have seen 
that at NLO C+ contains only on a single dimensionful 
scale t. In this appendix, we prove that this holds to all 
orders in perturbation theory. In Sec. IVII we argued that 
for a given parton channel {q — !■ qg, g — >■ gg^ or 5 — ^ qq) 



the matching coefficient C+ is a universal function, and 
our arguments here apply to all splitting channels. The 
main tool we use is RPI [20| . RPI restores the symmetry 
that is broken by choosing a particular basis of light-cone 
coordinates to label the collinear fields. For example, by 
choosing specific light-cone vectors n and n to define the 
2-jet operator O2, we have picked out a specific direction 
in an equivalence class of directions close to the jet di- 
rections that would all give the same results. RPI trans- 
formations correspond to reparameterizing to a different 
light-cone basis within the equivalence class, and invari- 
ance under such reparametrizations effectively restores 
the broken symmetry. In Ref. JJl] , RPI has been used in 
a similar way to what we follow here to understand the 
matching between operators with different multiplicities 
in SCET. 

There are three types of RPI transformations for each 
pair of basis vectors {n, n}. For our argument in this 
section we are only concerned with what is known as 
RPI-III, under which each direction n and its conjugate 
direction n transform as 

<^e"-<, nf^e— nf. (Al) 

Invariance under this transformation implies that any oc- 
currence of must be accompanied by an occurrence of 
. Each jet momentum q^ is RPI-III invariant, as it can 
be written as 



2K 



(A2) 



As discussed in Ref. [22l], the hard matching onto an A^- 
jet operator in SCET is RPI invariant. 
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CN{{uJi\){N\ON{{uJi\)\2). 



(A3) 

where uji = nt ■ qi, and we ignore the color structure as 
it is irrelevant for this discussion. Since the operator On 
is built out of gauge-invariant quark and gluon jet fields 
with given label momenta. 
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it transforms under RPI-III as {OjyduJi})) — ?> exp(ai -I- 
• • • -I- aN){ON{{^i})) ■ It follows that both the matching 
coefficient, CN{{uji}), and the operator matrix element 
together with the measure, (Ojv({<^i})) Yii ^i^ij sepa- 
rately RPI-III invariant. In what follows we suppress the 
integral over the uji labels. 

Consider the matching of On-i in SCET onto in 
SCET+, evaluated by calculating the 2 — >■ A^ matrix ele- 
ments of both operators with external momenta qi such 
that su/sij ^ Aj. In the SCET above the matching 
scale t — S12, the final states with momenta qi and 92 
are described by the same rit collinear sector with scaling 
Q(Aj , 1, Ai). Below this scale in SCET+, the final states 
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with qi and (72 are described by two separate cohinear 
sectors. The matching then takes the schematic form 
[see Eq. (|6J5l) ] 

{N\On-i\2)^C+{N\0+\N). (A5) 

where in pure dimensional regularization, C+ is given by 
the IR-finite terms in the matrix element (A^ | Oat- i| 2). 

After performing the BPS field redefinition, the matrix 
element can be factorized as 

{N\On^i\2) = (2|C„J0)(0|C„J1)---(1|C„„|0) 

x(0|yw-i|0), (A6) 

where Yiv_i stands for the product of A" — 1 usoft Wil- 
son lines. Since a separate RPI exists for each collinear 
sector, each collinear matrix element (IjCnjO) must sep- 
arately be RPI-III invariant. However, the only objects 
that each collinear matrix element can depend on are 

The only RPI-III invariant scale that we can define from 
these objects is q-f = 0. Therefore, the loop corrections to 
the collinear matrix elements are scaleless and hence zero 
in pure dim reg. A similar argument goes through for the 
usoft matrix element (0|Y/v_i|0). This is of course the 
familiar result that the virtual loop corrections to jet and 
soft functions vanish in pure dimensional regularization. 

If we apply a similar argument to the nt collinear ma- 
trix element (2|C„J0), we do not find that it vanishes in 
pure dimensional regularization. Since there are two par- 
ticles in the final state, the set of objects that the matrix 
element can depend on is 

{n-t, n'i, n'^, n^, n^, nt^, , q!^} . (A8) 

What RPI-III invariant, dimensionful parameters can be 
constructed out of these? While there are many choices, 
they fall into three categories: 

{91.2 = 0' 91-92, 91,2) 91,2)} , i = l,2,t. (A9) 

We now argue that all of the nonzero scales are the same 
order as t = 2qi • 92- To understand this scaling, note 
that all objects in Eq. (|A8|) can be constructed from the 
components of , q!^, and q^. Enforcing q^ = qf^ + q!^ 
and (Zi = ?2 = leaves only 6 independent variables in 
the components of qi, 92, and qt- We can choose these 6 
variables to be t, the energy Et = Ei + E2, the energy 
fraction z — Ei/ Et, where we are working in the regime 
where z is not near or 1, i.e., parametrically E1/E2 ~ 1. 
The last three variables are three angles. One is the 
azimuthal angle of the 1^-2 splitting, while the other 
two are trivial and merely specify the coordinate system. 
The dependence on the azimuthal angle is determined 
by the spin structure and therefore only appears as an 
overall dependence. 

In principle, there are two scales, t and Et, which the 
objects in Eq. (|A9I) can depend on, and t Ef. How- 
ever, using 

t = 2qi-q2 = {ni ■ qi)(ni • (72) = zEt{ni ■ 172) (AlO) 



it follows that fii ■ qi Et and ni • (72 ~ t/Et, and it is 
easy to see that all other nonzero • (71^2 and ■ qi 2 
have the same parametric scaling, rii ■ qi^2 t/Et and 
ni ■ qi^2 ~ Et- Therefore all possible RPI-III invariant 
combinations in Eq. (|A9p are of order t, while the Et 
dependence can only enter through higher-order power 
corrections of t/Ef. 

This scaling dependence then implies that at leading 
order in the power counting the only scale that all loop 
corrections can depend on is t, and after integrating over 
the azimuthal angle, the matching coefficient will only 
depend on t and z, as we have seen in the NLO result for 
C+ in Eq. (IS.lip . This is the desired result, valid at all 
orders in perturbation theory. 

Appendix B: All-Order Scale Sensitivity of the Soft 
Functions 

In Sec. I VII we derived the factorization theorem for 
the A^-jettiness cross section in the limit t ^ Sij. In 
this appendix we use RPI-III to derive the scale depen- 
dence of the csoft and usoft functions and show that 
they each depend on only a single parametric scale to 
all orders. The csoft function S+{ki,k2, n) only de- 
pends on the scale ^ ki/ \/st, and the usoft function 
S'Ar_i(fci, ^2, . . . , /cat, /i) depends on the scale Hsn-i 

We start with the csoft function, ^^(fci, k2, fi), defined 
in Eq. (j6.22p . and determine how it transforms under 
RPI-III. The structure of 5+ is independent of the par- 
tonic splitting channel, which only determines the color 
representation of the Wilson lines. The structure of the 
measurement function depends on the observable. Re- 
call that our discussion of the measurement in Sec. IIVBI 
applies to the geometric measure. 

In general, 5*+ depends on {ki,k2} through its mea- 
surement function in Eq. (j4.34l) and on the reference vec- 
tors {91, ^2, 9t}- The dependence on qi and q2 comes 
from the X Wilson lines and the measurement function, 
while the dependence on qt comes from the V Wilson line 
and the dependence on qt comes from the zero-bin sub- 
traction. The directions qi,2.t are in the same csoft RPI 
equivalence class and therefore transform under RPI-III 
as 

qi,2.t e"qi^2,t , Qi e.~°'qt . (Bl) 

To determine the transformation properties of 5*+ under 
RPI-III, we only need to consider the transformation of 
the measurement function since the Wilson lines are RPI- 
III invariant. The csoft measurement function transforms 
under the RPI-III transformation in Eq. (|Bip as 

M'%ki,k2) e-'" Y[ 5{h - 2q, ■ P^) (B2) 

i=l,2 

where we have rescaled the ki convolution variables to 
e^fci, which determines their transformation properties 
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under RPI-III. The analogous result also holds for the 
zero-bin measurement function. The csoft function there- 
fore transforms as S+ — ^ e~^"S'+. We can see this scaling 
in the RPI transformation of the tree level result 

^^^"(fci, fc2, /i) - S{ki) 5{k2) ^ e-2"J(fci) ^(fca) , (B3) 

which is simply the statement that the 5*+ has mass di- 
mension —2 and therefore has an overall scaling factor 
of l/(fci/c2) at all orders in perturbation theory multi- 
plying a RPTIII invariant function with nontrivial fci 2 
and ^ dependence. The convolution measure dfcidA:2 has 
a compensating e^" dependence and therefore the com- 
bination dfci dfc2 5'+ (A;i , /c2 , /i) is RPI-III invariant. It is 
useful to rescale the measure to factor out this 1/(^1/02) 
scaling from the csoft function, which effectively changes 
the measurement function to an RPI-III invariant mea- 
surement, 



x'=^(fci,fc2,M) = n '^(1 



1=1,2 



• Pf 



(B4) 




and similarly for the zero-bin measurement function. 
This allows us to write the csoft function in an RPI-III 
invariant form, while the combination dfci dfc2/(fci/c2)'S'+ 
is unchanged. We can now consider all RPI-III invari- 
ant scales that arise in the rescaled csoft function. This 
method of changing the convolution variable to write the 
delta function in an RPI-III invariant form is analogous 
to the method used in Ref. [l^ to analyze the RPI trans- 
formation properties in the hard matching onto SCET. 
The only RPI-III invariant dimension-one combinations 
that we can construct in the rescaled csoft function are 



(B5) 



where i,j S {l,2,i}, and recall that Sij = 2qi ■ qj. Addi- 
tionally, we can form the dimensionless RPI-III invariants 
Qt ■ cji and Sij/sik, but these are order 1, and so we can 
ignore them. Since all of the Sij are of order st = 2qi ■ q2, 
this shows that the only parametric scales that the csoft 
function depends on are k\^2/^fsti as advertised. 

The scale dependence of the usoft function Sn-i 
can be constrained using analogous arguments. The 
usoft function depends on the iV + 2 measure- 
ments fca, fcf,, fci, . . . , fc^r and the TV -I- 1 directions 
<?o , 9b , 9t , 93 , • • • (Zw ■ Note that all of these directions are 
well separated, so that any Sij — 2qi ■ qj ^ 1. Addi- 
tionally, in the ki and fc2 measurement functions the jet 
boundaries depend on the orientation of the nearby jets 
relative to the rest of the event, which is parameterized 
by an angle 0(. Since this angle is parametrically order 
1, we need not consider its dependence. 

The RPI-III transformation properties are more intri- 
cate now, as there are N + 1 separate coUinear sectors 
contributing usoft Wilson, each of which can in princi- 
ple have a separate RPI-III transformation parameter 
ai. The jet boundaries are not invariant under arbitrary 



RPI-III transformations, since the measurement contains 
comparisons of the form 



9{qi ■ pk < qj ■ Pk) , 



(B6) 



which do not have simple transformation properties un- 
der qi — > e^'Cji and qj — > G°''q_j- 

However, the usoft function has simple transformation 
properties under a restricted set of RPI-III transforma- 
tions with all ai set equal, so that all qi transform as 



e qt 



(B7) 



We call this a global RPI-III transformation, since it uni- 
versally affects all coUinear sectors. Note that in Sm~i 
there are no n directions to consider, since they do not 
enter in the Wilson lines or the measurement function. 
Under a global RPI-III transformation, the jet bound- 
aries are invariant, and the measurement function has a 
simple transformation property: 

MN{ki, . . . , fcAr) ^ e-(~+2)a_^^(^^^ . . . , few) , (B8) 

where ki e^ki under global RPI-III transformations, 
just as for 5*+. The usoft function therefore transforms 
as Sn~i — e~'^"'"^^"S'7v-i. This is compensated for 
by the transformation of the measure of the convolution 
variables under global RPI-III J], dfc, e(^+2)° JJ. dfe;, 
such that the combination Yii dfcj Sn-i is invariant. As 
before, it is useful to rescale the measure to factor out 
the Hi l/^i scahng dependence of Sn-i- The usoft mea- 
surement function then becomes 



Ml'i{k,})^'[[s(l 



■ pr 
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(B9) 

The rescaled usoft function is now invariant under 
global RPI-III. To determine its kinematic dependence, 
we construct all global RPI-III invariant dimension-one 
combinations, which are 
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for all z G a,b, 1,2, N and j,k S a,b,t, . . . , N with 
i ^ k. Unlike the S'+ function, there is no fii dependence 
in Sn-1- The only dimensionless RPI-III invariant com- 
binations of the qi that we can construct are the ratios 
Sij/ Ski, which are all order 1 and can again be ignored. 

We can further restrict the set in Eq. (jBlOl) using the 
fact that in the rescaled function, each ki only appears 
in the measurement function and always comes with qi 
in the combination (Ji/ki. (Since the qi can also arise 
from the Wilson lines, they can also appear without cor- 
responding factors of ki.) This requires either j = i or 
k = i in Eq. (|B10I) , reducing the set of possible combina- 
tions Sn~i can depend on to 
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TABLE II: Anomalous dimension coefficients for the various 
functions in tlie factorization theorem. 



Therefore, the only scales that the soft function can de- 
pend on to all orders in perturbation theory are those in 
Eq. (|Blip . Since all of the ^ sq , these scales are all 
of the same parametric order. 



Appendix C: Resummation Formulas 

In this appendix we collect all inputs needed for the 
resummation at NLL in Sec. lVIll The inputs required for 
the resummation at NNLL can be found for example in 
Appendix D.2 of Ref. [Hi and Appendix B.3 of Ref. [H]. 

The evolution kernels Ux for the different functions in 
Eq. (I7.3P which solve their evolution equations are writ- 
ten in terms of the evolution functions 

Vx{f^o,f^) ^ Cx'qr{f^o,fJ.) , (CI) 

where jx is the dimension of the evolution variable and 
Cx the coefficient of Fcusp (o^s) in the anomalous dimen- 
sion for each function, which are summarized in Ta- 
ble |lll The evolution functions Kr{nQ,fi), 77r(/^0:M), and 
^■yx if^o 7 A*) are given further below. 

The RG evolution of the hard function H2 following 
from Eqs. (15. 5p and (15. 6p is 

H2{Q\^l)^H2{Q^^lo)UH,{Q^^lo,^^), (C2) 



iO\''C2(Mo,p) 



Mo 



The RG evolution of the hard function H'^ following from 
Eq. (jSl^ is 



i?^ {t, z, n) ^ H1{t, z, ^0) {t, z, ^0, m) , 

-t — iO\''c_^(/io,M) ^ 
Mo 



(C3) 



Here, K^^{iiq, 11, x) explicitly depends on k and x via 

The solution for the soft function RGE in Eq. (|5.35p is 
given by [H 

= y dfc' S{k - k', pLo) Us{k', pio, A*) , (C4) 



Us{k,HQ,n) 
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Ks—7e VS 
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where Ks = Ksi^J,o,^J■), Vs = ?ys(Mo,M), and s = St or 
s — Sq for S+ or 5*2, respectively. The plus distribution 
is defined as 
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The RGE for the jet functions is given in Eq. (|5.18p . The 
solution has exactly the same structure as Eq. (IC4[) . 



(C6) 



J{s,fi) = J ds' J{s - s' , ^lo) Uj{s' , ^iQ, ^i) , 

gKj-lE rjj 







LMo Vo^ 





where Kj = Kj{fiQ, ^) and rjj = 77,7(^^0, m) as in Eq. (jCip . 
The convolutions of the plus distributions in the final 
cross section are easily evaluated analytically and can be 
found in Appendix B of Ref. [46| . 
The RGE functions are defined as 
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Expanding the beta function and anomalous dimensions 
in powers of as, 
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the explicit expressions of the evolution functions at NLL 
are 
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Here, r — as{^J.)/as{fJ,Q) and the running coupling is given 
by the three-loop expression 
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where X = 1 + Q;3(/io)/3o ln(Ai/Mo)/(27r). 

The coefficients of the beta function [54 
anomalous dimension [321 in MS are 
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The individual quark and gluon contributions to the non- 
cusp parts of the hard and jet anomalous dimensions in 
Eqs. (15.191) and (|5.43|) at one loop are 
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